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SUMMARY 



An extensive analysis was performed on 8.0 pb ^ of dimuon data produced in A 
collisions at 515 GeV/c to search for A}, events in the decay channel A|, — > J/\|/ A° , with 
J/\|/ — > 10,^ |i. and A° — > p 7U" (and for the conjugate reactions). The muon tracks from J/\|/ 
decays were refittedwith the mass constraint of two-body decay and with the constraint that 
both muon tracks intersect at a common point. The A*^ s were identified by their 
characteristic decay, giving the large fraction of their momenta to the protons. A° s were 
reconstructed in three regions of the E672/E706 spectrometer, in the SSD/target region, in 
the region between the SSDs and the dipole magnet, and inside the dipole magnet using an 
iterating algorithm. Several cuts were applied to the A° (and A^) to make a clean A^ 
(andA^) data sample. Among the cuts there was a K^^ mass cut, in which if the A*^ (or 
A°) had the mass of the K g under the hypothesis of both tracks being n 7U , the A 
(or A^) candidates were rejected. These give a clean A° (and A^) data sample of 575 ± 35 
A*^ (and A^) candidates. 



To search for the A^, — > J/\|/ A° (and charge conjugate reaction), J/\|/ s that passed 
the muon refit were combined with A*^ s (or A^ s) when they existed in the same event. 
The results show 2 events in the Aj, mass region. Using the E672 measurement of the bb 
cross-section, and considering the 2 A^, event candidates as signal with zero background, 

-2 

an upper limit to F(Ai,) * Br(Aj, J/V|/ A ) was found to be less than 6.2 xlO at 90 % 



xvu 



C.L. An upper limit was also calculate without using a K ^ mass cut for the A s (and 
A^s), and then, requiring that the J/\|/ s originate from secondary vertices, giving 
F(Ab) * Br(Ab ^ J/\|/ A° ) < 3.1 xlO"^ at 90 % C.L and F(Ab) * Br(Ab ^ J/\|/ A° ) < 3.2 
xlO"^ at 90 % C.L, respectively. 
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CHAPTER 1 



INTRODUCTION 

Since the discovery of the b quark in 1977 [1], much progress has been made in 
measuring and understanding the properties of beauty particles; this is primarily true for the 
By and B^j mesons, but not as much for the beauty baryons. The existence of (bud), the 
lightest of the baryons containing a b quark, has been somewhat controversial. In 1981, 
experiment R415, which used the Split Field Magnet (SFM) detector at the Intersecting 
Storage Ring (ISR) at CERN, reported an observation of the Aj, through the decay channel 
Ab ^ p D° 7l" [2] yielding a rest mass for the A^, of (5425!75^) MeV/c^. Later, in 1982, 
experiment R416, using an upgraded version of the R415 detector, performed the same 
search for A|, — > p 7t with a negative result [3]. This led to some discussion by the 
R415 collaboration claiming that experiment R416 did not have the necessary rejection 
power against charged hadrons; therefore, it should not have been able to observe beauty 
baryons at the ISR [4]. Further evidence for the Aj, was reported in 1991 by the former 
R415 collaboration in a second upgraded experiment, R422, at the ISR. Two different Aj, 
decay modes were observed: Aj, — > p D*^ 7t and — > Ag7t'''7t~7t~ [5]; the mass of the Aj, 
was found to be ( 5640^200 ) MeV/c^ and ( 5650!2^) MeV/c^ for each mode, respectively. In 
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1986, an experiment using the Fermilab Multiparticle Spectrometer reported the observation 
of a heavy baryon decaying into A!^K^'K^'K~Ti^n~ and having an invariant mass of 5750 
MeV/c^ [6]; however, this experiment did not claim to have observed the A^. CERN 
experiment UAl reported in 1991 the discovery of the decay channel Aj, — > J/V|/ A*^, 
claiming a signal of 16 ± 5 events above a background of 9 ± 1 events [7]; they measured 
the mass of the A^, to be (5640 ± 50 (stat) ± 30 (sys)) MeV/ c^. Since then, other 
experiments have searched for this decay mode such as ALEPH and OPAL at CERN [8], 
and CDF at Fermilab [9], however they have all failed to confirm UAl's observation. In 
1992, LEP experiments ALEPH, DELPHI, and OPAL reported evidence for the existence 
of the Aj, [10] through observation of the semi-leptonic decay channels Aj, — >A^1~ VjX 
and Ab A^ F Vj X in decays. 

In the fixed-target experiment E672/E706 at Fermilab, we have measured the bb 

total cross-section using our 1990 % -A interactions at 515 GeV/c [11]. In addition we 

+ ±0 0* 

observed exclusive decays of B hadrons such as B — > J/\|/ K and B — > J/\|/ K [11]. 

This suggests us to search for the A^, J/v|/ A^ decay mode, and to measure the product of 

the production fraction F(Ai,), times the branching ratio Br(Ai, — > J/\|/ A*^). 

The remainder of this chapter involves a review of the basic concepts in particle 
physics which are pertinent to the study of beauty baryons. It also presents detailed results 
of previous experimental searches for the A},. 



3 



1.1 Theoretical overview 



Present evidence indicates that matter is built from two types of fundamental 
particles called quarks and leptons, which are structureless and point-like on a scale of 
10"'^ m. Quarks carry fractional electric charge -e/3 or +2e/3, where e is the magnitude of 
the electron charge, and have spin 1/2. They come in several different flavors labeled u 
(up), d (down), s (strange), c (charm), b (bottom), and t (top). Each quark has it own 
internal quantum numbers. The u and d quarks are grouped in an isospin doublet with I = 
1/2, and with the third component I3 = +1/2 for the u and 13= -1/2 for the d. The s quark is 
assigned an internal quantum number called strangeness, with value S = -1. The c quark is 
assigned an internal quantum number called charm, with value C= -i-l. The b quark is 
assigned the bottom quantum number B= -1, and the t quark is assigned the top quantum 
number T = -i-l. Quarks also have a baryon number assigned to them, B = 1/3 for quarks 
and B = -1/3 for the antiquarks, which is conserved in any interaction. Because of then- 



Table I. The Quarks. 



Internal Quantum Electric Charge Rest Mass a 
Symbol _ 



Number 



: L doublet 



S = -l 
C=+l 

B =-1 
T = +l 



2-8 

5 -15 
100 - 300 
1000 - 1600 

4100-4500 
176,000; 199,000 



The quark masses are taken from Ref. [12], except the top 
quark mass which is taken from Ref. [13] and [14]. 
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masses, the u, d, and s quarks are referred to as light quarks and the c, b, and t as heavy 
quarks. Table I shows the electric charges and masses for the six quarks along with their 
internal quantum numbers. Each of the six quarks has its antiquark partner, which has the 
opposite quantum numbers and the same rest mass. 

The leptons carry integral charge, or ± e, and have spin 1/2. The neutral leptons 
are called neutrinos, and have very small (perhaps zero) rest masses. The leptons appear to 
come in doublets, with each neutrino being assigned a subscript corresponding to its 
charged partner. The three different types of charged leptons are known as the electron (e), 
the muon (/n), and the tau (t). Charged leptons are distinguished from antileptons by the 
sign of their charges. The neutrinos are longitudinally polarized: they have their spins 
opposite to their velocity vectors (left handed), while antineutrinos have spins in the same 
direction (right handed). A lepton number L^, L^, and of +1 is assigned to each type of 
lepton, respectively, and -1 to each t5^e of antilepton. The lepton number is always 
conserved in any interaction. The properties of the leptons are summarized in Table n. 



Table II. The Leptons 



Symbol 


Electric Charge 
Q/e 


Rest Mass ^ 
(MeV/c2) 


Antiparticle 


e~ 


-1 


0.511 




V 

e 





< 0.0051 




M~ 


-1 


105.6 




V 





<0.27 




i~ 


-1 


1777.1 


T + 


V 

T 





<31 





^ The mass values are obtained from Ref. [12]. 
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Quarks and leptons exist in three generations: The u, d, e ' and are the first 
generation; the s, c, fi , and are the second generation; and the b, t, t , and are the 
third generation. There are experimental measurements that indicate that there are only 
three generations \ 

Four fundamental forces (the gravitational, the electromagnetic, the weak, and the 
strong) govern the interactions between quarks and leptons. The gravitational force is by 
far the weakest of the four, and can be neglected for the study of interactions between 
elementary particles at tj^ical distances of the order of one femtometer. Thus, ignoring the 
gravitational force, the charged leptons have electromagnetic and weak interactions, while 
the neutrinos have only weak interactions. The quarks are subject to the electromagnetic, 
weak, and strong interactions. The field quanta, or mediators, for the electromagnetic, 
weak, and strong forces are the photon (7), the intermediate vector bosons and Z*^, and 



Table III. The Gauge Bosons 



Force 



Field Quantum 



Electric Charge Rest Mass^ 
Q/e (GeV/c2) 



Electromagnetic 



r 







< 3 X 10-27 



Weak 



W-, Zo 



±1,0 



80.22,91.18 



Strong 











^ The values for the boson masses are obtained from Ref.[12]. 
' See page 1333 of reference [12]. 
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the gluon (g), respectively. All of the field quanta have spin 1 and are called, in general, 
gauge bosons. Table III summarizes the properties of the gauge bosons. 

The electromagnetic force is described by the theory of Quantum Electrodynamics 
(QED). Here, the sources of the field are the electric charges of the interacting particles. 
The strength of the coupling constant of this force is called a and has a value of 
approximately 1/137. In the case of the strong force, which is described by the theory 
called Quantum Chromodynamics (QCD), quarks have an additional quantum number, 
whimsically called color charge, which is the source of this interaction. Each quark can 
have one of three primary color charges: red, green, or blue, denoted symbolically by R, 
G, and B, respectively. The antiquarks have, in similar fashion, color charges denoted by 
R,G,andB. Just as electric charge is conserved, so is color charge. Gluons also carry the 
strong-field color charge. They are bicolored objects and can be of eight different 
color/anti-color combinations^: 



Quarks combine to form colorless particles (color singlets) called hadrons, which are 
strongly interacting particles. In practice, we cannot observe isolated quarks or gluons; 
instead, we observe hadrons. The quark and gluon constituents of hadrons are generally 
known as partons. Hadrons that contain valence quarks are called mesons, and those 



RB , RG, BG , BR, GR, GB, 



RR-BB 



, and 



RR + BB- 2GG 




'With three colors and three anti-colors one expects 3^=9 combinations, but one of these is 
a color singlet and must be excluded. 
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that contain qqq (or qqq ) valence quarks are called baryons. The coupling constant of the 

2 

Strong interactions is called and it depends on the momentum transfer Q in a given 
reaction. Thus, is said to be a running coupling constant. 



1271 



(33-2nf)log 



qcd; 



where is the number of flavors, Aqqj is a parameter to be determined by experiment. 

2 2 2 

Therefore, if the /ly = 5 threshold is at Q = (2 mj,) =100 GeV , and with Aq^q = A5 = 
0.2 GeV, then is equal to 0.48. 

The weak force is described by the electroweak theory developed by Weinberg and 
Salam in the late sixties [15]. This gauge theory unifies the electromagnetic and weak 
interactions. The source of the weak force is a property called the weak charge, which 
gives quarks and leptons flavor-changing transitions. The corresponding effective 
coupling constant is called the Fermi constant, and has a value of Gp =1.2x10"^ GeV"^. 
The fundamental electroweak coupling constant, called g-^, is related to Gp through 
g^l%M-^ = Gp/ 42, where is the mass of the W boson. 
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The above brief introduction to quarks and leptons and the interactions among them 
is the modem view of fundamental particles and forces known as the Standard Model. 
Table IV illustrates Feynman diagrams^ of tj^ical interactions in the Standard Model. 



Table IV. Feynman diagrams of Typical Interactions. 



Interaction 



Charge 



Quarks 



Leptons 



Strong 



Color 




Electromagnetic Electric charge (e) 





Weak 



Weak charge(giy), 
giving u->d or Vg -> e' 
flavor-changing 



transions 





' For a description of how to read and construct Feynman diagram see Ref. [16]. 



1.1.1 Hadro production of b quarks 
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The hadronic production of b quarks is fundamental for the study of perturbative 
quantum chromodynamics (QCD) since results may be used to test and constrain the 
current form of the theory. The cross-sections for open-flavor (i.e., unbound) production 
of heavy quark pairs ( QQ ) are predicted by QCD. The lowest-order parton subprocesses 
are quark- antiquark fusion and gluon-gluon fusion, illustrated in Figure 1.1. 

The hadron-hadron (/^l/^2 QQX) production cross-section for heavy quark pairs 
is expressed as convolution of the parton-parton scattering cross-sections, G^j, with the 
distribution function of the partons, summed over all partons: 




Figure 1.1. Lowest-order QCD subprocess producing QQ heavy quark pairs. 
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c>(s,m,Q^) = Y^ I dx, ^dx2fi(xi,Q^)fj(x2,Q^)dij(S,m,Q^), (1.1) 

A /xiS 



where 5 is the square of the total center-of-mass energy in the hadron-hadron system, m is 

2 

the mass of the heavy quark, and Q is the mass factorization scale. The / and j indices run 
over all partons (quarks, antiquarks, and gluons) in the hadrons hj and respectively. 
The variable Xj is the fractional four-momentum vector of the parton in hadron hj, and X2 
is that of the parton in hadron Thus, in the parton-parton system the square of the 
center-of-mass energy is ^ ~ ^jXjS. / (-^i^G^) dxj is equal to the probability that parton i 
in hadron carries a fraction of S, analogously, fj{^X2,Q^) dx2 for parton j in hadron 
h2. These parton densities are empirically determined from electron and neutrino deep- 
inelastic scattering (see Ref. [17]), and the parton-parton scattering cross-section, 
(7, ^ (s,m, (2^), can be calculated in perturbative QCD [18]. 

The partonic cross-sections in a next-to-leading-order (NLO) or 0{al^ power- 
series expansion are usually written in terms of so-called scaling ratios as 



where ri = — ^-1 or ri = ^^. The dimensionless functions ^ ^ (?/) and fijiXl) 
4m s 

represent the Bom approximation and the Of a^j corrections, respectively, and ^j^(?7) 
appears when the mass factorization scale, , deviates from the square of the heavy flavor 




(1.2) 



4m^ 
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quark mass, m^. NLO calculations of the partonic cross-sections have been computed by 
Nason, Dawson, and ElUs (NDE) [19]. The bb total cross-section for hadron-hadron 
collisions at current fixed-target energies has been calculated to NLO using NDE results by 
Berger [20], and also by Mangano, Nason, and Rodolfi (MNR) [21]. Experimentally, the 
bb total cross-section has been measured at fixed-target energies by experiments WA78 
[22] and NAIO [23] at CERN and E653 [24] and E672/E706 [11] at Fermilab; Figure 1.2 
shows the results for 7t~N collisions and also the theoretical predictions of Berger and 

MNR for several different mass factorization scales, Q^, and heavy-flavor quark masses, 
m^. It is important to mention here that in hadron-hadron collisions at collider energies 

such as that in CDF and DO (4s = 1800 GeV) the bb total cross-section is three orders of 
magnitude larger than at our fixed-targert energies [25]. 
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Figure 1.2. The bb total cross-section for 7l~N collisions as measured by references. [22], 
[23], [24], and [11], and the theoretical prediction by Berger [20] and MNR [21]. The 
uncertainties shown for the measured values are the statistical and systematic contributions 
added in quadrature. 



1.1.2 Fragmentation of b quarks 
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After production of QQ pairs, fragmentation, or hadronization, occurs subsequently. 
During this process the color forces organize the QQ pairs with other quarks created in this 
process into colorless hadrons. Fragmentation is governed by soft, non-perturbative 
processes that cannot be calculated from first principles and can only be modeled. In the 
LUND string-breaking model [26] for example, the color energy stored in the color field 
increases as the individual QQ quarks separate; then, at some point, the color field has 
enough energy to produce a (light quark pair) or a (light diquark pair) and the 

string breaks. This effect is shown in Figure 1.3. In this model the ratio of the probability 
for gggg production to that for is 0.09 : 1.0, where the probability for qq production 
approximately 0.3. [27]. Ifab (heavy) quark combines with a single light quark, the result 
is a beauty meson; if it combines with a light quark pair, the result is beauty baryon. The 




Heavy Baryon 



Q Heavy Antibaryon 

Figure 1.3. Schematic picture of fragmentation in which breaking of the lines of the color 
force between the separating heavy quarks of a QQ pair produce a light llll diquark pair 
which combines with the heavy quarks to form a heavy baryon pair . 
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probability F(Ab) of a b quark fragmenting into a A^, or into a beauty baryon that 
subsequently decays into a Aj,, has been measured to be 7.2 ± 1.2 % [28]. 

1.1.3 b-hadr on decay 

Several kinds of Feynman diagrams contribute to the non-leptonic decays of b 
hadrons; examples are shown in Figure 1.4. The processes shown are called spectator 
(external or internal), exchange, annihilation, and penguin (electromagnetic or hadronic) 
diagrams. In spectator diagrams the b quark decays in to a c quark (shown), or u quark 
(not shown) and a W~ boson which, in turn, decays primarily into ud or cs (favored in 
the Cabbibo-Kobayashi-Maskawa (CKM) scheme), or into e~y^, v^, or T~v^. In these 
decays, the lighter quarks in the b hadron act only as spectators and do not participate in the 
decay. It is important to mention the difference between "external" and "internal" spectator 
diagrams. In the external diagram, the colorless W~ is allowed to decay into three possible 
ud or cs color singlet final states. In contrast, in the internal diagram, in order to form 
colorless final states only one of the three possible virtual W~ decays is allowed. 
Therefore, the internal diagram is suppressed by a factor of 1/ with respect to the 
exterior diagram, where A^^ is the number of colors. 




Electromagnetic Penguin Hadronic Penguin 



Figure 1.4. Feynman diagrams for the non-leptonic decay of a b hadron. 
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The total b-quark decay width from the spectator diagrams, F,^"'^^, has been 
computed to be [29] 



.totSP _ <J T^h 



192%' 



2.72|VbJ' + 6.92|VbJ' 



(1.3) 



where G is the Fermi coupling constant, m^ the mass of the b quark, V^c and V^u are the 
CKM matrix elements that represent the — > c and b^u quark transitions, respectively. 



In the W-exchange and annihilation decays, both initial-state quarks are involved in 
the weak vertex; the classic example is 7C~ — > |J.~v^ (or tt" — > e~Vg). For b-hadrons these 
types of decays are heUcity-suppressed. To zeroth order in QCD, the total hadronic width 



for W exchange, T^^^, is given in Ref. [29] to be 



rb^-(bq^qiq3) = 



3 87U 



'qqs 



fbqmb(mi +m|)M(mim3,mb), ( 1.4) 



and that for annihilation is 



Fb— (bq^q2q3) = 3 



91 

87t 



V, 



bq 



V, 



qiqa 



fbqinb(m^+m^)M(m2,m3,mb), (1.5) 



where the Vqg' are CKM matrix elements. 



M(mi,m3,mb) = 



21 2 , 2\ 

mb(mi +m3j 



1- 



mj + m3 
V ^b J 



1- 



m^-mj 
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and f^q, called the decay constant, has dimensions of mass and is related to the amplitude 
of the bq wave function at the origin by f|,q =12|\|/(0)| /M^q (where Mi,q is the mass of 
the bq system. The value of |v|/(0)| can be calculated in a similar fashion as for the 
hydrogen atom. It involves the reduced mass of the system, m(,mq /(m(, + mqj ~ mq, and 

|2 ^ ^ 

the coupUng constant a^; that is, |\|/(0)| oca^mq. Thus, by taking Mj,q~mi, there 
follows 



-iWex+annih { ^ \^ 



1-totalSP "'s 



mq 



(1.6) 



Since widths are proportional to rates, this indicates that the ratio of the (W exchange + 
Annihilation) decay rate to the total spectator decay rate is, in general small, since m^, » 



mq. 



The decay width for the penguin diagram is given by [29] 



rP-(bq^sq) = ^(c5 + 3c6f |-|Vt,f|%3ff2qm3, ( 1.7) 



where Cj and Cg are constants with values of 0.02 and -0.04, respectively, defined in Ref. 
[29]. The penguin contribution to b-hadron decay is negligible (a few per cent) with 
respect to the spectator contribution. 
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In the spectator-type decay where b ^ ccs, the. cc can bind to form a r]^ , J/v|/, Xi, 
¥' ,or a higher-mass charmonium state. The relative ratios for the production of some of 
these states is predicted in Ref. [30] to be 0.57 : 1.0 : 0.27 : 0.31 for t]^ :J/\|/: Xi:¥' . 

The — > J/\|/ A° decay has no nonspectator contributions; thus, the theoretical 
computation of the branching ratio for this decay is relatively clean. The theoretical 

-3 

prediction is that the Br(A^ J/\|/ A ) is less than 10 [3 1]. 

1.2 Evidence for the existence of the 

While b-flavored baryons have long been predicted by the quark model, only very 
recently (1992) has their existence been confirmed. Evidence for A^, baryons was 
reported by LEP experiments ALEPH, DELPHI, and OPAL [10] through observation of 
the semi-leptonic decay channels Aj, — > A^ 1~ Vj Xand Ay, — > A"^ 1~ Vj X in Z" — > bb 
events. The analysis was based on correlations between the baryon and lepton produced in 
the decays, since b baryons are expected to produce A^ 1 , A^l"*" or A^ 1~ , A~ 1"*" pairs, 
and not A^ 1^, A^l~ or A^l", A^l^ pairs (see Figure 1.5). The experiments found an 
excess of events of the correct baryon-lepton combination pairs, and interpreted this excess 
as A), decays. The results, in the form of the product of production fraction and branching 
ratio are, summarized in Table V. 
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Table V. Results on A -lepton and Aj.-lepton correlations in Z decays at LEP. The first 
uncertaiaty is statistical and the second is systematic. 





F(A,,) *Br(Ai, ^ aVvjX) 


F(Ab)*Br(Ab^A+rViX) 




(10"^) 


(10'') 


ALEPH 


6.1 ±0.6 ± 1.0 


1.51 ±0.29 ±0.23 


DELPHI 


3.0 ± 0.6 ± 0.4 


1.18 ±0.26 


OPAL 


2.91 ±0.23 ±0.25 


0.83 ± 0.28 



1.3 Previous sear ches for the J/\|/ A decay channel 

The quark diagram for the — > J/\|/ A° decay channel is shown in Figure 1.6. 
This diagram is an internal spectator diagram. Thus, it is colored-suppressed, since the cc 
pair form a color-singlet state in order to be bound. Only CERN experiment UAl has 
observed A^, s in the J/\|/ A^ decay mode. They found a signal of 16 ± 5 events above a 
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background of 9 ± 1 events [7], and measured the mass of the Aj, to be (5640 ± 50 (stat) ± 
30 (sys)) MeV/c^ and F(Ab) * Br(Ai, J/\|/ A^ ) to be (1.8 ± 0.6 (stat) ± 0.9 (sys))x 10"l 
ALEPH, OPAL, and CDF have failed to confirm UAl's signal (see Refs. [8] and [9]). 
The results found in their search are summarized in Table Table VI. 




Table VI. Results from other searches for the A(, — > J/\|/ A^ decay channel. 
Experiment F( A^, ) * Br( A^ ^ J / \|/ A° ) 

UAl (1.8 ± 0.6(stat) ± 0.9 (sys)) x 10-3 

CDF < 0.5 X 10-3 

ALEPH < 0.4 X 10-3 

OPAL < 3.4 X 10-4 
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1.4 Analysis objective 

The object of this thesis is to search for the A^, beauty baryons in the A^, — > J/\|/ A^ 
channel (and the charge conjugate reaction) using the 1990 data from fix-target experiment 
E672 at Fermilab. The J/\|/ s are identified by their decay into |i. and the A^ through 
the decay channel A^ — > p 7t . By relating the number of A^, — > J/\|/ A^ events to the E672 
measurement of the b-quark cross-section obtained from the same data, the production rate 
F(A^,) * Br(A|, J/\|/ A^ ) will be extracted. (Throughout this thesis, whenever a state is 
mentioned its charge conjugate state is also implicitly implied.) 



CHAPTER 2 



THE MWEST SPECTROMETER 

The MWEST spectrometer, was located at the meson-west area at Fermilab. It was 
used simultaneously by experiments E672 and E706. The experiments accumulated data in 
1990 and 1991. The output signals from all detector systems were available to both 
experiments. E672 concentrated on dimuon events, whereas E706 was interested in events 
that contained direct photons [32]. Experiment E672 was designed to study the 
hadroproduction of high-mass muon pairs which resulted from heavy-quark production, 
particularly charmonium [33]. This chapter begins with an overview of the complete 
spectrometer and then presents a detailed description of the individual components of the 
apparatus that were used to obtained the data for the analysis in this thesis. 

2.1 Overview of the detector 

The physical layout of the combined E672/E706 apparatus is shown in Figure 2.1. 
A right-handed orthogonal coordinate system was associated with the spectrometer, with its 
origin near the target and with the z-axis aligned with the beam direction. The x and y axes 
were in the horizontal and vertical directions, respectively. The apparatus began with a 
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hadron shield used to absorb beam halo and a set of scintillation-counter veto walls used to 
veto events with halo particles at the pretrigger level. Downstream was a Cu and Be target, 
a tracking system composed of a set of 16 silicon-strip detectors (SSDs), a dipole magnet 
with average p^ kick of 0.4457 GeV/c, 16 proportional-wire-chambers (PWCs), and 4 
straw drift-tube detector planes. Following this was a liquid-argon electromagnetic 
(EMLAC) and hadronic (HADLAC) calorimeter, and a forward calorimeter for energy 
measurement of forward-going electrons, photons, and hadrons. Finally, about 20 meters 
downstream of the target was the muon spectrometer consisting of a concrete hadron shield 
(beam dump) to prevent hadrons produced in the underlying event into the muon system, a 
toroid magnet with average p^ kick of 1.3 GeV/c, 12 muon PWCs, and two scintillator- 
hodoscope planes (muon hodoscopes HI and H2) with associated electronics for 
generating a trigger signal upon the detection of two muons. Each of the elements pertinent 
to the analysis in this thesis are discussed in more detail. 



Beam 
Hodoscope 

Veto Walls 

Target/ 
SSDs 



LAC 



beam 



PWCs 



I Neutron l_i_J V 
1 f Hadron Absorber a ^^^^^ 



\ \ 
y J 



Z Shield 



Dipole 
Magnet 



Tubes 



+ 



Muon PWCs 
|jA, jiB |J.1-|J.4 




Beam 
Dump Toroid 
Magnet 



Muon Hodoscopes 
HI and H2 



m 



10 m 



20 m 



30 m 



Figure 2. 1. The physical layout of the MWEST spectrometer. 
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2.2 Beam 

The Fermilab Tevatron operated on a 57.2 second cycle. During which, 23 
seconds were used for the spill of the beam to the experimental areas. The Tevatron 
provided an 800 GeV/c proton beam to the switchyard, the area where the beam was split 
and directed to the three major beam areas: meson, proton, and neutrino. The beam had a 
19.7 ns bucket structure from the characteristic radio frequency of the Tevatron. 

2.2.1 Beamline 

The beam used in the experiment was produced in the following way: The portion 
of the beam sent to the meson west area, which had a t5^ical intensity of about 2 x 10^^ 
protons per spill, collided with a beryllium target of dimensions, 46.5 cm in length and 
2.22 cm in diameter, corresponding to 1.14 interaction lengths. This produced a secondary 
K~ beam of 515 ± 3 GeV/c average momentum [34], and intensity at its maximum of 
2x10^ K~ per spill [35]. The same bucket structure applied to the secondary beam, where 
the probability of a bucket being occupied by a single particle was about 10%, and the 
probability of being occupied by more than one particle was approximately 2%. 

Both pion and proton beams were delivered to the meson-west area through the 
Fermilab meson-west beamline, which was built specifically for the experiment. This 
above ground beamline was designed to transport positively and negatively charged 
particles with a momentum of up to 1 TeV/c. 
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2.2.2 Cerenkov counter 

The secondary beam used by the experiment in 1990 was not a pure pion beam, 
97.0% of the beam particles were 7r~s, 2.9% ^"s, and 0.1% antiprotons [35]. A 
Cerenkov detector was used to tag the beam particles. The detector had dimensions of 
42. 1 m long and radius of 24.4 cm, and it was positioned in the direction of the beam. 
Helium gas, at pressure between 4-8 psi, was used as the radiator. For a more detailed 
explanation of the Cerenkov counters see reference [35]. 

2.3 Hadron shield and veto walls 

A large stack of iron, of dimensions 4.7 m long, 4.3 m wide, and 3.7 m high, with 
a hole in its center (for the beam particles to pass through), was placed at the end of the 
beamline, just inside the experimental hall. This served as a hadron shield to absorb the 
beam halo, those particles traveling in the beam direction, but not on the beam axis. Also 
"spoiler" magnets in the upstream beamline were employed to sweep halo particles away 
from the beam. 

Halo muons, produced primarily from 7t~ decays from the beam, could however, 
pass through the iron shield very easily, and subsequently throughout the entire 
spectrometer. To tag the muons so that they would not start the dimuon trigger (for 
example, one muon from the target and a halo muon), a series of three "veto walls" 
scintillation counters were used. One was placed just upstream of the hadron shield (VW3) 
and the other two just down stream (VWl and VW2), see Figure 2.1. The veto wall VW3 
consisted of 18 rectangular scintillator counters, and it covered an area of about 3 m x 3 m. 
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Veto walls VWl and VW2 each consisted of matrix of thirty two 50 cm x 50 cm scintillator 
counters, also covering an area of approximately 3 m x 3 m. These two walls had a 10 cm 
offset in the x-y plane to cover the gaps from one another. The three walls had a hole in the 
central region to allow the beam to pass without generating a signal. A coincidence 
(VW1+VW2) * VW3 between the walls was established to veto events containing halo 
muons, and this was done at the pretrigger level. 

2.4 Target 

The experiment used a target composed of two pieces of Cu and two pieces of Be, 
as shown in Figure 2.2. The two Cu targets were 0.8 mm thick, and were separated by a 
0.26 cm air-gap. Downstream of this, 0.55 cm, was the first piece of Be, which had a 
thickness of 3.71 cm, and was followed by a 1.02 cm air-gap. Finally, there was the 
second piece of Be, which was 1 . 12 cm thick. 

2.5 Upstream tracking system 

This tracking system, located near the target, was used to reconstruct the tracks 
(trajectories), and the vertices made by the tracks of the outgoing charged particles accepted 
by the spectrometer. It had four major parts: the SSDs, the upstream PWCs, the dipole 
magnet, and the STRAWS. The STRAWS will not be discussed here, since they were not 
used for this analysis. 
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2.5.1 Silicon-strip detector system 

The arrangement of the SSD system is shown in Figure 2.2. It consisted of 16 
separate planes, 6 of them upstream of the target and 10 planes downstream. The SSD 




Figure 2.2. The physical layout of the target and SSD system. 
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planes measured the position of a charged particle by putting out a signal induced by the 
particle on the strips etched onto the siUcon wafers. The 16 planes were setup in pairs, 
consisting of one plane with its silicon-strips oriented in the x-direction (horizontal) and the 
other in the y-direction (vertical). The 3 x-y planes upstream of the target, had an area of 3 
cm X 3 cm, and a pitch of 50 |i.m. The 5 pairs of planes downstream of the target all had a 
5 cm X 5 cm area, and except for the plane just after the target, which had a pitch of 25 |i.m, 
they all had a 50 |J,m pitch. 

The vertex position resolution of the SSD system was of 10 |xm in the x- and y- 
direction, and 350 mm in the z-direction. The angular resolution was 0.1 milliradians. 
This tracking device played an important role in the reconstruction of secondary vertices, 
those produced from decays in flight of particles emerging from the primary interaction. 

2.5.2 Dipole magnet 

A conventional liquid helium dipole magnet was used to measure the momentum of 
the charged particles. This was done by measuring the deflection of their trajectory as they 
passed through the magnet. The dipole magnet was operated at a current of 1050 A, giving 
rise to an approximately uniform magnetic field of magnitude 0.61 15 Tesla. The magnetic 
field was oriented in the negative y-direction, and it extended through 2.416 m (length of 
the magnetic field). A mirror plate was placed in the front (upstream side) and back of the 
magnet to reduce the fringe field. Each plate had a rectangular opening to allow particles to 
pass through, the one in front was 35.56 cm wide x 25.4 cm high, and the one in the back 
127.0 cm X 91.44 cm. A polyethylene bag filled with helium gas was installed in the 
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central region of the magnet to reduce the multiple scattering of the practices when 
traversing through the magnet. 

Since the direction of the magnetic field was aligned with the y-direction (vertical), 
the trajectories of the charged particles were only deflected in the horizontal x-view (x-z 
plane). The effective change, in the x-component of the particle's momentum was by the 
so called " p-p ^^^.^^ of the magnet, 

PTkick = 0.3BoL, (2.1) 

where Bg, the magnitude of the field, is in Tesla, and L, the length of the field, is in 
meters. This was equal to 0.4457 GeV/c in the x-direction. 

2.5.3 Upstream proportional wire chambers 

Downstream of the dipole were a set of 16 proportional wire chambers (PWCs). 
They were grouped into four modules, named: PWCl, PWC2, PWCS, and PWC4. Each 
module contained a plane with its anode wires aligned with the x- (horizontal), y- (vertical), 
u- (rotated + 37° from the vertical), and v- (rotated -53° from the vertical) directions. The 
physical dimensions of the four modules were different in order to obtain a constant 
angular acceptance. PWCl was 1.63 m wide x 1.22 m high, PWC2 and PWC3 were 2.03 

2 2 

m , and PWC4 was 2.44 m . The anode signal planes were made of gold-plated 0.8 mm 
in diameter wires, spaced at 2.54 mm. For each anode plane there were two cathode 
planes, one on each side, separated by a distance of 5.74 mm. These cathode planes were 
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made of graphite-coated 1.0 mm thick mylar sheets. The chambers were filled with a gas 
mixture of 18% isobutane, 2.2% isopropyl alcohol, 0.1% freon, and 79.7% argon. The 
PWC system is described in grater detail in [36] 

2.6 Muon system 

The muon system was located at the down stream end of the spectrometer about 20 
meters from the target. The physical layout of the muon detector is shown in Figure 2.3. 
It consisted of: 8 upstream muon PWCs, a beam dump (made of steel and tungsten), a 
conventional toroid magnet, four downstream muon PWCs, two 16-segment scintillating 
counter hodoscopes, and concrete walls to shield the muon PWCs. 

2.6.1 Upstream muon proportional wire chambers 

The first two muon PWCs named |i.A and |iB were installed to reduce the number 
of triggers caused by the halo muons and to improve the linking tracks between the 
upstream tracking system and the muon system. Both )J,A and |lB consisted of four anode 
signals planes: x, y, u, and v. The x and y planes were oriented in the horizontal and 
vertical directions, respectively. The u planes were rotated 10° from the vertical, and the v 
planes were rotated -10*^ from the vertical. The anode planes were made of 1 .0 millimeter 
diameter gold-plated wires, separated 3.05 mm. The dimensions of the active areas and 
the numbers of wires for each anode plane are listed in Table VII. 
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Muon hodoscopes HI and H2: 
16-segment scintillator- 
hodoscope planes 



Downstream muon PWCs: 




Upstream muon PWCs; Each with 4 anode planes, 
horizontal, vertical and ±100 from vertical of 0.3 cm pitch 



Figure 2.3. The physical layout of the muon spectrometer 
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The cathode planes consisted of 1650 3.5 mil diameter Cu/Be wires, separated 1.02 mm. 
The distance between adjacent cathode and anode planes was 0.95 cm. Two ground 
planes were also included in each PWC module, one in front of the first cathode plane and 
the other after the last cathode plane. The ground planes were made of 1725 3.5 mil 
diameter Cu/Be wires, separated 1.02 mm [37]. The gas mixture used in the PWCs was 
76 % argon, 15 % isobutane, 8.9 % methylal, and 0.1 % freon. For more information on 
these PWCs see reference [37]. 



Table VII. Upstream muon PWC specifications. 


PWC plane 


Orientation from 
the vertical 


Number of wires 


Active area 
(m2) 


|J.A X 


900 


256 


1.7 X 0.8 


|XA y 


00 


560 


1.7x1.7 


|lA u 


IQO 


352 


1.7x1.1 


|i,A V 


-IQO 


352 


1.7x1.1 


[lB X 


900 


544 


1.7 X 1.7 


\lB y 


00 


544 


1.7 X 1.7 


|i,B u 


100 


448 


1.7 X 1.3 


^lB V 


-100 


448 


1.7x1.3 



2.6.2 Toroid 



An iron polarized toroid magnet was used to provide a second measure of the 
momentum of the muons. The toroid was 2.44 m long, had an outer radius of 1.35 m, an 
inner radius of 16.8 cm at the upstream end, and an inner radius of 19.7 cm at the 



33 



downstream end. The magnetic field created by the toroid (operated normally at 1700 
amperes) was in the azrmuthal direction. It had a magnitude of 2.24 T at the inner radius 
and 1.74 at the outer radius. The field had no measurable dependence in the azrmuthal 
angle. Muons received a momentum kick in the radial direction which decreases linearly 
from 1.64 GeV/c at the inner radius, to 1.27 GeV/c at the outer radius. The measurement 
of the momentum of the muons was used at the trigger level and was later compared to 
momentum measurement of the upstream tracking system. 

2.6.3 Muon pr etrigger hodoscopes 

Two scintillating hodoscopes were used in the muon spectrometer, arranged as 
shown in Figure 2.3. The muon hodoscopes both had a radius of 1.5 m to match the active 
area of the downstream muon PWCs. They consisted of 16 triangular scintillator segments 
ordered in a circular petal pattern, as shown in Figure 2.4. Each of the segments was made 
of PS- 10 plastic scintillator encased in aluminum sheets. The scintillating light produced 
upon the passage of a muon through a scintillator segment was collected by a BBQ doped 
waveshifting bar and re-emitted at a wavelength of 410 nm. This light was collected by 
RCA 8575 photomultipUer tubes operating at 2100 volts. At this voltage the hodoscopes 
were at least 90 % efficient. The daisy counters were used at the pretrigger level to tag the 
muons produced in the underlying event . 
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2.6.4 Downstr earn muon pr oportional wir e chambers 

The main components of the muon system were the four downstream PWCs, 
named: |J,1, [l2, \x3, and [l4. These chambers detected the trajectories of the muons as they 
passed throughout the muon system. Their arrangement is shown in Figure 2.3. The 
physical dimensions of PWC |il are shown in Figure 2.5 (a), and the specifications of all 
the downstream muon PWCs are listed in Table VIII. Each chamber contained two 
cathode and one anode signal planes (see Figure 2.5(b)). The cathode planes were made 
of copper clad G-10 epoxy-fiberglass sheets of 1.6 mm thickness, etched to give a pattern 
of parallel strips 9.5 mm wide and separated 1 mm. The strips were oriented in the u- 
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(+ 45 from the vertical) and v- (- 45 from the vertical) directions (see Figure 2.5(c)). The 
anode planes were oriented in the y-direction (vertical) and made of 25 |im diameter gold- 
plated tungsten wires, separated by 6 mm, except for |i,2 which had 4 mm pitch. The 
anode wires of |il, |i3, and |J,4 were connected in groups of three adjacent wires to each 
amplifier card, and in pairs for |i.2, so the effective pitch of each PWC was 1.2 cm. The 
detail explanation of the muon PWC readout, and in general the data acquisition system 
(DA), is explained in the next chapter. 



Table VIII. Downstream muon PWC specifications. 
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Anode-Cathode separation 


all 


9.5 mm 
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Figure 2.5. Downstream muon proportional wire chamber: (a) physical dimensions; (b) 
cross section of the chamber; (c) cutaway view of the chamber. 



CHAPTER 3 



DATA ACQUISITION SYSTEM 

The data acquisition system of the MWEST spectrometer recorded the events of 
interest for both experiments, E672 and E706. E672 triggered on events that contained at 
least one high-mass dimuon, and E706 had several LAC based triggers [38]. This chapter 
describes the data acquisition system of the whole detector, however, only the muons 
system readout is described in detail. 

3.1 Overview of the data acquisition system 

The data acquisition system (DA) was divided into four parts, three CAMAC based 
PDP-11 systems and one FASTBUS system. A schematic diagram of the DA is shown in 
Figure 3.1. The muon system was readout using a PDP-11/34, named MU. It 
communicated with the muon system CAMAC interface through a Jorway 411, which read 
the CAMAC crates on a serial CAMAC highway. The upstream PWCs, SSDs, and 
CERENKOV counters were readout through a PDP-11/34, named NEU, and the forward 
calorimeter (FCAL) through a PDP-11/34, named ROCH. The FASTBUS system was 
used to readout the LAC and straw drift-tubes planes (STRAWS). The PDP-lls and 
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FASTBUS were connected to a device developed by Fermilab, called a Bison Box [39], 
which controlled the beginning and ending of the spill interrupts, as well as the event 
interrupts (triggers). The PDP-lls were also connected to a DEC Micro VAX 11 computer, 
using communication device (CD) links [40]. The Micro VAX was running the 
VAXONLINE software system [41], which concatenated the readout data of the PDP-1 Is. 
This data, along with the FASTBUS data, constituted a complete event. THE VAXOLINE 
also controlled the beginning and ending of the data taking process, kept a dynamical event 
pool on the Micro VAX, for on-line data monitoring, and wrote the collected data to 8 mm 
exabj^e tapes 

3.2 Muon system readout 

The electronic readout of the muon system consisted of several CAMAC crates 
modules which included a LeCroy PCOS HI DataBus Interface (LeCroy 4299). This 
collected data from the muon chambers, ADCs, TDCs, scalars, the dimuon trigger 
processor (DTP), and sent it to the PDP-1 1 MU (see Figure 3.1). 

3.2.1 Muon pretrigger hodoscope readout 

The output signals from the 16 phototubes of each hodoscope were sent through a 
RG58-C coaxial cable to two 16-channel CAMAC multiplicity logic units [42]. The two 
modules discriminated the received signals and set a latch for each channel that exceeded a 
certain preset voltage threshold during an external set gate. The multiplicity of the tracks in 
the muon hodoscopes was then calculated from the multiplicity of the latch bits. The 
CAMAC module could be preset to a certain multiplicity threshold, so if the calculated 
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multiplicity was grater than or equal to the set threshold, a NIM level was output. During 
normal data taking the multiplicity threshold was set to two for both muon hodoscopes. 
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Figure 3.1. Diagram of the MWEST detector data acquisition system. 
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3.2.2 Upstream muon proportional wire chamber readout 

All the muon chambers were readout using a LeCroy PCOS HI system 
(Proportional Chamber Operating System) [43]. The readout signals were sent to LeCroy 
2731 latch modules, which resided in CAMAC crates. Each crate could hold 23 latch 
modules, and one controller (LeCroy 2738). The controller module, which queried each 
latch module, stored the addresses of the active channels in an internal buffer, then sent 
them to the data acquisition system. The controller also had an ECL port which transferred 
data ( at a rate of 10 times that of CAMAC) to the trigger logic circuit. 

3.2.3 Dimuon trigger processor 

The dimuon trigger processor (DTP) was used to do an on-line calculation of the 
dimuon invariant-mass. It consisted of 7 CAMAC double-width modules. These received 
the data output of the muon PWCs PCOS III. The data was sent by the LeCroy 2738 crate 
controller through the fast ECL port. A flow chart diagram of the DTP is shown in Figure 
3.2. The wire-hit data information from the chambers, |J,1 and |l4, was sent to two 
modules named POINT. These modules computed in parallel (using the wire hits) the 
space points of the trajectories of the muon tracks. The location of these spatial points was 
checked against a look-up table stored in Programmable Read Only Memory (PROM) to 
assure that the points were in the feducial volume of the muon spectrometer. The list of 
valid points was stored in a Random Access Memory (RAM) and sent to two other 
modules named TRACK. 
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The TRACK modules also received the wire hit information from the |l2 and |i3 
muon PWCs. Using this additional information and the valid space points, they calculated 
in parallel the possible associated tracks in the x (x-z plane) and y (y-z plane) views. This 
list of tracks was sent to two other modules named MOMENTUM and TARGET. The 
TARGET module, which also received the hit information from the upstream muon 
chambers [lA, and |i.B, checked to see if the tracks were consistent with the hits found in 
the |jA and |J,B. This was done in the following way : An imaginary straight line was 
constructed from the center of the target to the point in space were the track in question met 
the center of the toroid; the imaginary line was compared with the wire-hits of |i,A and |a,B, 
and if the hits on at least 3 out of 4 planes in each chamber were consistent with the 
imaginary line, the associated track was said to be a valid track. Again, in parallel with the 
TRACK module, the MOMENTUM module calculated the momenta of the tracks from 
their bend in the toroid assuming they originated from the target. The last step of the DTP 
was to compute the invariant-mass of every pair-combination of valid tracks. This was 
done (assuming the tracks were from muons) by a module named MASS, which received 
the list of valid tracks and their associated momenta . If at least one pair-combination of 
valid tracks gave an invariant-mass greater than a preset threshold, a success signal was 

2 

output. In the 1990 run, the mass threshold was usually set between 0.5 GeV/c and 1.0 
GeV/c^ 
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3.3 Dimuon trigger 

The E672 dimuon trigger was activated after the detection of both a beam particle 
and an interaction, and if the DA was not busy. It was done in two stages: Level 1, or 
pretrigger; and level 2. 

The pretrigger required that the hits in the muon hodoscopes be consistent with a 

track multiplicity in the muon detector of least two tracks. If this was true and no particles 

were detected by the vetowalls, level 2 was started. Muons produced in the target required 

approximately 15 GeV to penetrate the spectrometer material and reach the muon 

hodoscopes, all particles reaching them were assumed to be muons. The average hit 

.4 

multiplicity was 2.3. The pretigger rate was 2 x 10 per interaction. The pretigger 
efficiency for two muons penetrating the system was 0.76 and remained constant over the 
data taking period. 

Level 2 was based on the requirement that at least one of the dimuons had invariant- 
mass above a certain mass threshold. The invariant-mass calculation was done by the DTP 

based on the information in the muon chambers. If at least one muon pairs satisfied the 

2 

DTP conditions the event was recorded to tape. The threshold of 0.7 GeV/c resulted in a 
trigger rate of 2 x 10 ^ per interaction.. The average DTP processing time was 10 |xs per 
pretigger, which included 5 |J,s to decode the muon chamber data. The combined efficiency 
of the chambers and the DTP algorithm was 0.77 for dimuon events. 
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Figure 3.2. Schematic diagram of the dimuon trigger processor. 



CHAPTER 4 



EVENT SELECTION AND RECONSTRUCTION 

The raw data were written onto 8 mm exabyte tapes. These contain the complete 
information from every piece of apparatus in the entire spectrometer for all the different 
events that where triggered (dimuon and LAC triggers). This chapter describes the data 
reduction from the raw data to the data sample used in the analysis of this thesis. 

4.1 Event selection 

The first step towards reducing the data was to extract the dimuon triggers from the 
raw data tapes. This step reduced the amount of the data to about 18 % of the complete 
data sample, yielding approximately 5 million dimuon triggers (corresponding to a 
luminosity of 8 pb V To further reduce the data, only information from the upstream 
tracking system and the muon system were extracted. The information from the LAC and 
forward calorimeter was excluded. By excluding this portion of the data, the event size 
was reduced by a factor of 8. The set of events that were extracted is known as the "MNS" 
sample. 
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4.1.1 Dimuon preselection 

Because of inefficiencies of the veto walls, some halo muons activated the dimuon 
trigger. Also, muons produced by beam interactions with the forward calorimeter could set 
off the dimuon trigger. This happened often since the forward calorimeter only had a 3.2 
cm diameter hole in the center. If it was not properly aligned, the beam could possibly 
interact with the calorimeter and produce muons. To avoid events that might have any of 
the two problems just mentioned, the MNS dimuon sample was run through a preselection 
program. The program required that at least two muons with a momenta greater than 20 
GeV/c each be present in the event. Both muons had to be consistent with originating from 
the target region. Also, the track segments from the muon system had to be linked with an 
upstream PWC track segment. These requirements reduced the MNS sample by 86 % . 
The set of events that survived the requirements of the dimuon preselection program is 
known as the "DIM" sample. This sample consists of about 750 thousand events. To 
reduce the data even more, the events were required to have at least one opposite sign 

2 

dimuon pair with an invariant-mass greater than 2.0 GeV/c . Approximately 35,000 events 
survived these requirements. This set of events is called the "PSI" sample. It is the actual 
data set used in the analysis of this thesis. 

4.2 Event reconstruction 

The event reconstruction program was written by Prof. Jack Martin [44]. The 
program reconstructed the "tracks" (trajectories) of the charged particles detected by the 
spectrometer. It first found the track segments detected by the SSD, upstream PWC, and 
muon PWC systems. Using the track segments, it linked the tracks from the various 
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systems: The SSDs with the upstream PWCs, these are labeled SSD-PWC linked tracks; 
and the upstream PWCs with the muon PWCs, these are called PWC-MUON Unked 
tracks. It also linked the track segments throughout the complete detector, these tj^e of 
tracks are called SSD-PWC-MUON linked tracks or fully linked tracks. Finally, with the 
reconstructed tracks it searched for the location of the vertex in the event. All these steps 
are described in detail below. 

4.2.1 Track finding 

The track finding algorithm for the muon PWC system began by searching each of 
the view planes: u, v, and y, for groups of adjacent wire hits. These groups of hits were 
called hit clusters. The center of each hit cluster (with its corresponding uncertainty of half 
the width of the cluster) was used as a view hit. Having the list of view hits in each view, 
a search for single-view tracks was performed. They were found by doing a search of all 
the possible roads containing only one hit cluster per plane in at least 3 of the 4 planes. A 
road was defined as a quadrangular portion of the active area of a given plane view of the 
muon PWC system. It extended from the first PWC (|il) to the last PWC (|i4). Figure 
4. 1 shows the first road used in the algorithm, and the subroads into which it is divided. 
Roads that contained more than the allowed number of hit clusters were further subdivided 
using a scheme as in Figure 4.1. The algorithm was repeated until the roads satisfied the 
conditions mentioned above. The roads that did not meet the criteria were disregarded. 
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As a simple example of the track finding algorithm, consider the hit cluster pattern 
found in the muon system as shown in Figure 4.2. After the first iteration, it is easy to see 
that only subroad 3 has no more than one hit per plane in at least 3 out of the 4 planes. In 
this case, it has hit clusters on all the planes. This road is not subdivided further since it 



First road 




Original road 




Subroad 1 



Subroad 2 



Subroad 3 



Subroad 4 



Figure 4.1. Illustration showing the first roads used in the track algorithm and the scheme to 
subdivide a road into other roads. 
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now meets the track finding criteria. All the other roads are disregarded. Of course the hit 
cluster pattern found in a real event is far more complex than this simple example. 

Once the final set of roads is found, a straight-line fit is performed to the hit clusters 
within each road. This fit determines the slope and intercept of the tracks in each view. 
After this, a three-dimensional space-track finding is performed in a similar fashion as in 
the single-view track finding. Space-roads are constructed and checked to see if any 
single-view tracks are consistent with belonging to the space-road. The space-roads that 
contain more than one view-track are further subdivided. This is done until a set of space- 
roads which only contain one single-view track per view is found. The hit clusters 
belonging to each space-road are then fit to a straight line. The fit gives the parameters of 
the space-track as the slopes and intercepts of the track in the x-z and y-z planes (x and y 
views, respectively). 

The track finding for the upstream PWC and SSD systems are essentially the same 
as for the muon system. The only difference is that when fitting the hit clusters to find the 
final parameters of the tracks, single wire hits are used instead of the hit clusters (one wire 
hit from each cluster). The parameters of the single wire hit track combination, that give 

2 

the smallest X from the straight-line fit are kept as the parameters of the track. This is 
done because the high multipUcity of the wire hits gives rise to spuriously large hit cluster 
widths. No space-track finding is done for the SSD system, since it only had two plane 
views. 
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Figure 4.2. Illustration of a simple example of the track finding algorithm, the arrow 
shows the muon track. 
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4.2.2 Track linking 

After reconstructing the track segments in each of the tracking devices, these are 
linked from one device to another. The track linking fit determines the momenta of the 
tracks by using the bend in the dipole magnet, or in the toroid, for the case of muon tracks. 

The SSD-PWC track linking is performed in the following way: Track segments 
from the SSD and upstream PWC are extrapolated to the center of the magnet in each view. 
The segments from the y-view must meet one another within the error of track segments. 
The segments from the x-view must meet each other within 5 mm. The linking in this case 
is done in each view, separately, because the SSD system only has single-view tracks and 
no space-tracks. The reason for the linking requirements to be different in each view is 
because the dipole magnet only deflected the tracks in the x-view. Thus, the y-view tracks 
should have approximately the same parameters as seen in the SSD and upstream PWC 
systems. Figure 4.3 illustrates the SSD-PWC linking. Track segments that meet the 
linking criteria are fit using the slopes and intercepts of each track segment: a, b, e, f; and c, 
d, g, h (see Figure 4.3). The y-view fit is done with the constraint that the slope and 
intercept of the track are the same in both the SSD and upstream PWC systems. In the x- 
view, the fit outputs the momentum of the linked track (asl/P), the slope, and the intercept 

2 

of the track as seen by the SSD system in the x-view. Tracks that pass the % quality cut of 
the fits are said to be SSD-PWC linked tracks. PWC space-tracks can have up to three 

2 

SSD-x and three SSD-y link combinations if they pass the X cut. These ambiguities are 
resolved at a later stage. 
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Figure 4.3. SSD-PWC linking scheme. 
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The PWC-MUON track linking is done in the following way: The track segments 
from the upstream PWC and muon PWC systems are extrapolated to the center of the 
toroid. Track segments that meet one another with in 12 cm in both the x and y views are 
considered as linked track candidates. Figure 4.3 shows the PWC-MUON linking scheme. 

2 

The linked track candidates are fitted using a standard % minimization algorithm, using the 
parameters of the track segments: e, f, g, h, i, j, k, and 1. Also, the momentum (in the 
form 1/P) and two multiple scattering angles and ^y) are included in the fit. Candidates 

2 

that pass a % quality cut are said to be PWC-MUON linked tracks. Again, in the case of 
ambiguities, up to 3 PWC-MUON linked combinations are kept for each muon track. 
Matching of these tracks to the SSD-PWC linked tracks and to the vertex position are used 
to resolve the ambiguities at a later stage. A global fit is done to the muon tracks. These 
are linked throughout the entire detector. To do this, the SSD-PWC linked tracks are 
extrapolated to the center of the toroid and matched with the muon tracks in the same 
manner as the upstream PWC segments, described above. Here, the fit is performed using 
all the parameters of each segment composing each track: a, b, c, d, e, f, g, h, i, j, k, and 1. 
Also, the momentum of the track in the dipole ( written as 1/P), the momentum in the toroid 
( written as 1/P), and two multiple-scattering-angles (j)^ and (j)y are included (see Figure 

2 2 

4.3). Once again, a % minimization algorithm is used. Track combinations that pass a % 
quality cut are called SSD-PWC-MUON linked tracks, or fully linked muon tracks. 

To improve the resolution of the SSD-PWC linked tracks, the SSD segments 
associated with a vertex were refitted using the vertex position as an additional pseudo- 

2 

plane hit. A % fit like the one for the PWC-MUON linked tracks, is now performed on the 
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SSD-PWC tracks, yielding the momenta and the track parameters in each view at the 
production vertex. 

The fractional uncertainty in the momentum measurement of a fully linked (SSD- 
PWC-MUON Unked) muon track is: 



^ = 0.0005 p; 
P 



as seen in Figure 4.4. The over all efficiency to reconstruct a fully linked muon track is 
86%. 
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Figure 4.4. The fractional uncertainty of momentum vs. momentum for fully linked muon 
tracks. 
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4.2.3 Vertex finding 

The vertex finding algorithm begins by looking at the SSD x and y views for 
intersections of two track segments. In each view, intersections that occur approximately at 
the same position (within errors) are combined. A list of possible multi-segment 
intersections is made for each view. The two lists are compared, and only intersections that 
have consistent z positions (within errors) in both views are kept. This results in a three 
dimensional position (x, y, z) of the vertices and a list of the tracks associated with each 
one. 

A fit minimizing the impact-parameter of each track connected to its vertex is 
performed to improve the resolution of the measurement of the vertex position. Figure 4.5 
shows the distribution of the primary vertices found in the events of the dimuon PSI 
sample (descried in section 4.1.1) for those dimuon events that are tagged as J/\|/s\ This 
distribution clearly shows the target configuration as well as the SSD planes (see sections 
2.4 and 2.5.1 for the target and SSD plane configurations). The comparison with the 
actual physical size of the target indicates a resolution of approximately 3 mm in the 
longitudinal direction. Figure 4.6 shows the transverse impact-parameter of the SSD-PWC 
linked tracks to the primary vertex for the events in the PSI sample. This distribution 
demonstrates a resolution of about 10 |xm for the measurement of the vertex position in the 
transverse plane (x-y plane). 



' See chapter 5 for the definition of J/\|/ events. 




Figure 4.5. Distribution of the primary vertices found in the "PSI" sample for those 
dimuon events that are tagged as J/\|/ s. 



56 




5 10 15 20 25 30 35 40 45 50 



Transverse impactparameter (|a.m) 



Figure 4.6. Transverse impact-parameter of the SSD-PWC linked tracks with the primary 
vertex, for events in the "PSI" sample. 
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4.3 Characteristics of reconstructed dimuon events 

E672/E706 has studied the hadroproduction of charmoniuom by measuring the 
\|/(2S), %j, X2' directly produced J/V|/ fractions contributing to the inclusive J/\|/ 
production (see Ref. [45] and [46]), and compared these values with those predicted by the 
theoretical charmoniuom models: color evaporation [47] and color singlet [48]. We have 
also, as mentioned previously, measured the bb cross section by using J/\|/ s from 
secondary vertices [11]. This measurement will be used later to extract the production rate 
for Ab s, F(Ab)*Br(Ab J/\|/ A*^). Another E672/E706 study, in which the author has 
been directly involved is the measurement of the atomic number (A) dependence of the 
production cross sections for vector mesons p/co, (|), and J/\|/ [49]. The results of this study 
are not included in detail here in order to make this thesis coherent. The following is a brief 
summary of the results of these three studies. 

The invariant mass spectrum of our 1990 opposite sign muon pairs is shown in 
Figure 4.7. This is our so called "DIM" sample (see 4.1.1 for the for its definition). This 
sample contains approximately 36 K p/to, 10 K (j), and 12 K J/\|/. A closer look at the J/\|/ 
region (see Figure 4.8) shows a \|/(2S) signal which contains 270 ± 35 (stat) ± 50 (sys), 
\|/(2S) ^ |i events. The \|/(2S) is observed in both the |i and J/v|/ n modes. 
Figure 4.9 shows the J/\|/ n invariant-mass spectrum. A fit to the signal yields 224 ± 
44 (stat) ± 20 (sys), V|/(2S) J/\|/ %' events. Differential cross sections for J/\|/ and 
\|/(2S) can be found in [45]. The total integrated cross sections were measured to be 
Br( J/\|/ -^ii'if) a(7t~ Be J/\|/ + X)/A = (9.3 ± lA(sys)) nb/nucleon for Xp > 0.1, and 
Br( \\f(2S) J/\|/ %^ if ) o(7t" Be \\f(2S) + X)/A = (7.4 ± 1.5 (stat) ± 1.2 (sys)) 
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nb/nucleon for ;cp > 0.1. The fraction of inclusive J/V|/ yield due to V|/(2S) meson decays 
was found to be 0.083 ± 0.017 (stat) ± 0.013 (sys) [45]. 




Figure 4.7. Invariant-mass spectrum of the opposite sign muon pairs. 
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Figure 4.8. Invariant- mass spectrum of the opposite sign muon pairs in the J/\|/ region. 
The solid line is a fit to the data. 
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Figure 4.9. Invariant-mass spectrum of J/\|/ % % combinations. The solid line is a fit to 
the data. 
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Xi and X2 states are reconstructed via their decay into J/\|/ y. The y s are detected 
via their conversions into e^ e pairs. Figure 4.10 shows the J/\|/ y spectrum for converted 
y s. Clear peaks are seen for both Xi and X2- A fit to the signals yields 57 ± 13 {stat) ±16 
{sys) and 40 ± 10 {stat) ± 14 {sys) Xi and X2' respectively. The integrated total cross 
sections associated with these values are Br( Xi J/^ y) <s{n Be ^ Xi + X)/A = (50.5 ± 
9.0 {stat) ±4.1 {sys)) nb/nucleon for Xp > 0.1 and Br( X2 ^ J/¥ Y) ^(^ Be ^ X2 + X)/A 
= ( 35.5 ± 7.8 {stat) ± 7.6 {sys)) nb/nucleon for jcp > 0.1. Thus the fraction of inclusive 
J/\|/ yield due to Xi and X2nieson decays are 0.24 ± 0.04 {stat) ± 0.03 {sys) and 0.017 ± 
0.04 {stat) ± 0.02 {sys), respectively. 

Combining the 0.083 ± 0.017 {stat) ± 0.013 {sys) of J/\|/ s due to \|/(2S) and 
assuming 0.02 of the inclusive J/\|/ are from Xo and b-hadron decays, the fraction of the 
total inclusive J/V|/ produced directly is 0.49 ± 0.05 {stat) ± 0.05 {sys) [46]. A comparison 
of our measurements with the theory are shown is Figure 4.1 1. 




Figure 4. 10. Invariant-mass spectrum of J/\|/ y combinations for converted y s. The solid 
line is a fit to the data 
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Figure 4.11. (a) A comparison of our measurements of the fraction of inclusive J/\|/ yield 
due to \|/(2S), %i, %2' directly produced J/\|/ s, to the fractions predicted by the color 
evaporation model [47]; and (b) to the fractions predicted by the color singlet model [48]. 
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E672 has also measured the bb cross section [11]. This analysis was performed by using 
J/\|/ s emerging from secondary vertices in the regions of the target where only air was 
present. E672 found an excess of 8 events in the air-gap regions of the target, and 
attributed these to b ^ J/Vj/ + X decays. Which correspond to an inclusive cross section of 
75 ± 31 istai) ± 26 {sys) nb/nulceon for all Xp. Some exclusive B hadron events were also 
reconstructed. Figure 4.12 shows the combined invariant-mass spectrum for B* — > J/\|/ 

±0 0* 

K and B J/\|/ K events. There are five events near the nominal mass of the B, 3 of 

± 0* 
them are from J/\|/ K and 2 from J/\|/ K . 

In another study, we look at the atomic number dependence of the total and 
differential cross sections for the vector mesons p/co, (j) and J/\|/ and low-mass Drell-Yan 
pairs in the Feynman-x region 0.1 < Xp < 0.8. By parameterizing the total cross-section as 
Go A*^, we find the values for a to be equal 0.74 ± 0.01 {stat) ± 0.02 {sys), 0.80 ± 0.01 
{stat) ± 0.02 {sys), and 0.92 ± 0.02 {stat) ± 0.02 {sys), for these vector mesons, 
respectively. We find the value of a for the Drell-Yan dimuon continuum equal to 1.16 ± 
0.08 {stat) ± 0.02 {sys) for dimuon masses between 4.0 GeV and 7.0 GeV. We found no 
significant dependence of the a parameter on xp or the transverse momentum p^ [49]. 
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Figure 4.12. Combined invariant-mass spectrum for B — > J/\|/ K and B — > J/\|/ K 
(hatched). 



CHAPTER 5 



THE J/\|/ ^ |i SAMPLE 

The J/V|/ [i^ \x event selection begins by including all dimuons from the "PSI" 
data sample (the procedure to reduce the raw data to the "PSI" sample is explained in 

2 

section 4.1.1.). Each dimuon is refitted using a standard % technique, and only dimuons 

2 

that survive the % cut are kept, but no more than one dimuon per event. This chapter 
describes the requirements for the selection of J/\|/ \i events used in the Aj, — > J/\|/ 
analysis of this thesis. The details of the dimuon refit are explained in Appendix A. 

5.1 Selection of the J/\|/ \i' sample 

The dimuon invariant-mass distribution for the "PSI" sample is shown in Figure 
5.1. A fit to this distribution that included the Monte-Carlo resolution functions for the J/\|/ 
and V|/(2S), and two exponential functions for the background, yields a J/\|/ mass of (3.097 

2 

± 0.001) GeV/c . This value of the J/\|/ mass is in agreement with the world average value 
from the particle data group [12]. The J/\|/ signal region is defined as the mass interval 

2 2 

from 2.85 GeV/c to 3.35 GeV/c . This region contains about 13,053 muon pairs, with a 
background subtracted signal of 1 1,500 J/\|/ events. 
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The muon tracks in the signal region are then refitted with the constraint that both 

2 

muons intersect at a common pomt and that the invariant-mass is equal to 3.097 GeV/c . 
This constrained fitting technique improves the momentum resolution as compared to initial 
track finding by a factor of 2. Also, the vertex position resolution of dimuons improved by 
15 %. See appendix A for plots of the momentum and vertex residuals before and after the 

2 

fit. The quality of the fit was used to reduce combinatorial background. A cut on the % 
per degree of freedom of the fit, to be less than 5.0, reduced the |i background in the 
J/\|/ signal region by 50%. A total of 12,340 muon pairs survived the refit criteria. 
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Figure 5.1. The dimuon invariant-mass distribution for the "PSI" sample. The solid line is 
a fit to the data that includes the Monte-Carlo resolution functions for the J/\|/ and V|/(2S), 
and two exponential functions for the background. 



CHAPTER 6 



THE A" ^ p 71 (AND CHARGE CONJUGATE) SAMPLE 

The great majority of s produced in this experiment decayed after the 
target/SSDs. Therefore, the tracks of their daughter decay products (proton and pion) were 
not detected before their trajectories were deflected by the magnetic field of the dipole 
magnet. Thus, their reconstruction is more complex. The — > p 7t decays can be 
reconstructed if the A*^ s decayed upstream of the PWC system, see Figure 2.1. In 
particular there are three regions in the spectrometer where this reconstruction can be 
achieved: near the target/SSDs, upstream of the dipole magnet, and inside the dipole 
magnet, see Figure 6.1. In the region near the target/SSDs, the proton and pion are 
detected by the SSD and PWC systems, but, for the other two regions they are only 
detected by the PWC system. This is important, because the standard way of determining 
the momentum of a charged particle is by measuring the deflection of its trajectory upon 
passing through a magnetic fleld. This chapter describes the procedure and algorithms 
used to reconstruct A^ — > p 7t decays. Also, to confirm the validity of the reconstruction 
algorithms, k'^^ — > tc^ tc decays (which because of their lifetime have the same topology 
as A^ — > p 7U decays) are reconstructed in similar fashion, and they are also used to clean 
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the A*^ (and A*^) sample. I remind the reader that the throughout this thesis whenever a 



particle or a reaction is mentioned the charge conjugate particle or reaction is implicitly 
implied. 
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Figure 6.1. Regions in the spectrometer where A s can be reconstructed. 
6.1 Reconstruction algorithms for the A *^ 



An important characteristic of A decay in the laboratory frame, is the imbalance in 
momentum between the decay proton and pion. Because of the large difference in mass 
between the the proton and pion, the proton from the A^ decay will have on the average 
larger momentum than the pion'. In the a'^ ^ p 7U events accepted by the E672/E706 



' Mp = 0.938 GeV/c^ , 0.1389 GoWlc and in the center-of-mass frame p* = 0.101 

GeV/c, thus, Ep= 0.943 GeV and Ejj= 0.172 GeV. Therefore, after boosting the proton 
and pion momenta to the laboratory frame, on the average, the proton wiU have a larger 
momentum than the pion by a factor of approximately Ep / {- 5.48). 
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detector, the proton always had larger momentum than the pion (this was determined by 
Monte-Carlo). Thus, when performing the two-track combinations to reconstruct the A°, 
the track with the largest momentum of the two is always assumed to be the proton. 

6.1.1 Reconstruction near the target/SSDs region 

The topology of the A° — > p 7U decay in this region is illustrated in Figure 6.2. For 
these decays, the tracks of the daughter particles are reconstructed using the SSD and PWC 
systems. Their momenta are measured by linking the SSD tracks to their proper PWC 
tracks as explained in section 4.2.2. The A° s are then reconstructed by taking two SSD- 
PWC linked track combinations, assigning the proton mass to the track with the largest 
momentum of the two, and the pion mass to the other one. The proton and pion tracks are 
required to form a vertex. This means that the distance of closest approach between the 
two tracks is less than 50 |im. In order to see a clear A° signal, this vertex should be at 
least 1.5 cm downstream of the production vertex of the A*^, which in this case is assumed 
to be the J/\|/ decay vertex. According to Monte-Carlo, an impact parameter of a track to the 
primary vertex of less than 200 |i.m means the track is associated with the primary vertex. 
— At least 1.5 cm 1 



J/V|/ decay vertex 5 / A° decay ^^^if^ .0 ^ 

assumed to be the ^ ^ ^ ^ ^^^^ ^ , — — ■ A reconstructed 




p track 
a'' recoi 

a" production vertex ^ '-•^^^^^^^'^ P n trajectory 



TT track 

Figure 6.2. The illustration shows the topology of a A^ ^ p 7U decaying near the 
target/SSD region. 
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Therefore, the impact parameter of the pion 5^ relative to the A production vertex is 
required to be greater than 200 |im. Since the proton has a tendency of following the 
direction of its parent A^, it is easy for the vertex reconstruction program to assign the 
proton track to the primary vertex. The cut on the impact parameter of the proton 5p relative 
to the A° production vertex is therefore looser than that for the pion and it is required to be 
greater than 100 |xm. Also, according to Monte-Carlo the impact parameter 5^o of the 
reconstructed A*^ relative to its production vertex should be less than 120 |im. Using all 
these criteria a clear signal for A^ (andA^) is seen in Figure 6.3. The A^(andA^) signal 

2 2 

region is defined as the mass interval from 1.112 GeV/c to 1.118 GeV/c , and there are 
10 A'^(and A'^) candidates. A fit to this distribution that included a gaussian function, and a 
second order polynomial function for the background, yields 9 ± 3 (stat) background- 
subtracted events, and a A°(andA^) mass of (1.114 ± 0.0004) GeV/c^, with a FWHM 
mass resolution of 2.3 MeV/c^. This value of the A*^(and A*^) mass is in agreement with the 
world average value from the particle data group [12]. 
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Figure 6.3. Invariant-mass distribution using only SSD-PWC tracks for p % (and p 7U^) 
track combinations intercepting in the target/SSDs region. The solid line is a fit to the data. 
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6.1.2 Reconstruction in the region upstream of the dipole 
magnet 

Taking advantage of the two-body decay kinematics, and assuming the production 
vertex of the s to be the J/v|/ decay vertex, the A*^ s that decayed in this region can be 
reconstructed. The first step is to determine the A^ decay vertex (x, y, z). The dipole 
magnet only bent the trajectory of the electricly charged particles in the x-z plane. 
Therefore, the y and z coordinates are found from the intersection of the proton and pion 
PWC tracks in the y-z plane, see Figure 6.4. Knowing y and z, and assuming an x 
coordinate for the A° decay vertex, the corresponding momenta of the proton and pion can 
be computed. This can be determined by extrapolating the PWC tracks in the x-z plane to 



y-z plane 



target/SSD (y,z) 
□ 




Dipole magnet 
x-z plane 



target/SSD (x,z) , ^ ' 



Dipole magnet 



PWCs/STRAWS 




76.93 cm 197.73 cm z 
Figure 6.4. Illustration of the proton and pions PWC tracks, for those A^ s that decayed in 
the region upstream of the dipole magnet. The solid lines are the reconstructed PWC 
tracks. 
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the center of the dipole magnet, and using the A decay vertex to estimate the deflection of 
their trajectories throughout the dipole magnet, see Figure 6.4. See also Appendix B for 
the explicit calculation of the momentum vectors. 

To find the actual x coordinate of the A^ decay vertex, the assumed position of the x 
coordinate is scanned across the x-axis for the value such that the vector sum of the 
transverse momentum of the proton and pion with respect to the direction of flight of the A*^ 
divided by the sum of their magnitudes (this ratio is called the " relative p^ ", see Figure 
6.5) is a minimum. Real A^ decays should have, in principle, minimum relative p^ equal to 
zero; thus, cutting on this quantity removes undesired background to A*^ decays originating 
from accidental intersections of tracks. Figure 6.6 shows the distribution of minimum 
relative p-j- for proton and pion tracks intersecting upstream of the dipole magnet. 
According to Monte-Carlo studies, 99 % of reconstructed A*^ s have minimum relative pj 
less than 0.4. Thus, a cut on minimum relative pj of less than 0.4 is imposed on the data. 

A'' direction 




Figure 6.5. The illustration shows the topology of a A — > p 7C decay, and the definition 
of relative p-p. 



76 




50 - 



25 - 



0.1 0.2 0.3 0.4 0.5 0.6 

minimum relpt 

Figure 6.6. Distribution of minimum relative pp for tracks intersecting upstream of the 
dipole magnet. 



To reject more undesired background to A , the proton and pion PWC tracks are 
required to not be Unked to any SSD track. Also, because of the small difference between 
the A^ and proton mass, the proton tends to follow the trajectory of the parent A^. Thus, 
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the impact parameter of the proton 5p relative to the A production vertex is required to be 
less than 1.5 cm (this value was determined from Monte-Carlo). The invariant-mass 
distribution for the p 71 (and p 7t^) combinations for this region is shown in Figure 6.7. A 

clear signal for A'^(and A*^) is seen. The A*^ (and A*^) signal region is defined as the mass 

2 2 

interval from 1.105 GeV/c to 1.125 GeV/c , and there are 479 candidates. A fit to this 
distribution that included a gaussian function, and a fourth order polynomial function for 
the background, yields 220 ± 23 (stat) background subtracted events, and a A*^(andA'^) 
mass of (1.116 ± 0.0004) GeV/c^, with a FWHM mass resolution of 8.0 MeV/c^. This 
value of the A^(and A^) mass is in agreement with the world average value from the particle 
data group [12]. 
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Figure 6.7. Invariant-mass distribution using only PWC tracks for p 7U (and p 7U^) 
combinations for tracks intercepting in the region upstream of the dipole magnet. The solid 
line is a fit to the data. 
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6.1.3 Reconstruction inside the dipole magnet region 

For those A^s that decayed inside the dipole magnet the reconstruction algorithm is 
similar as for those that decay upstream of the magnet. Again, the reconstruction algorithm 
begins by finding the (y, z) coordinates of the A^ decay vertex. The y and z coordinates 
are determined by the intersection of the proton and pion PWC tracks in the y-z plane, see 
Figure 6.8. The determination of the x coordinate is the same as before. Knowing the y 
and z coordinates and with a given x coordinated for the A^ decay vertex, the momenta of 
the proton and pion are estimated by measuring the radii of their unique circular trajectories 
through the dipole magnet projected onto the x-z plane. These circular trajectories are such 
that they include the point (x, z). Also at the downstream end of the dipole magnet (z= 



y-z plane 
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197.73 cm 318.53 cm 

Figure 6.8. Illustration of the proton and pion PWC tracks, for those A*^ s that decayed in 
the region inside the dipole magnet. The solid lines are the reconstructed PWC tracks. 
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318.53 cm) they match the slopes of the reconstructed PWC tracks (see Figure 6.8). See 
also Appendix C for the detailed calculation of the momentum vectors. 

To find the actual x coordinate of the A*^ decay vertex, the assumed position of the x 
coordinate is scanned across the x-axis for the value of x such that the relative p^ is a 
minimum (see Figure 6.5 for the defenition of relative p^). Real A° decays should have, in 
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Figure 6.9. Distribution of minimum relative pj for tracks intersecting inside the dipole 
magnet. 
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principle, minimum relative pj equal to zero; thus, cutting on this quantity removes 
undesired background to decays originating from accidental intersections of tracks. 
Figure 6.9 shows the distribution of minimum relative p^ for proton and pion tracks 
intersecting inside the dipole magnet. Again, according to Monte-Carlo studies, 99 % of 
reconstructed A°s (andA^s) have minimum relative p^ less than 0.4. Thus, a cut on the 
minimum relative p^ to be less than 0.4 is imposed on the data. 

To further reduce the background to A^, the proton and pion PWC tracks are 
required to not be linked to any SSD track. The impact parameter of the proton 5p relative 
to the A^ production vertex is required to be less than 1.5 cm. The invariant-mass 
distribution for the proton-pion combinations that intersect this region is shown in Figure 
6.10. A clear signal for A°(andA°) is seen. The A^ (andA^) signal region is defined as 

2 2 

the mass interval from 1.10 GeV/c to 1.13 GeV/c , and there are 582 candidates. A fit to 
this distribution that included a gaussian function, and a fourth order polynomial function 
for the background, yields 346 ± 26 (stat) background subtracted events, and a A°(and A^) 
mass of (1.115 ± 0.0006) GeV/c^ with a FWHM mass resolution of 14 MeV/c^. This 
value of the A^ (andA^) mass is in agreement with the world average value from the 
particle data group [12]. 
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Figure 6.10. Invariant-mass distribution using only PWC tracks for p 7t (and p 7U ) 
combinations for tracks intercepting in the region inside the dipole magnet. The solid line 
is a fit to the data. 
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6.2 The k", ^ n signal 

Since the meson has a lifetime of the same order of magnitude (10 s) as the 
A°, the — > 71^ 71 decays have the same topology as the A° — > p 7i decays. Thus, the 
^ 7U decay becomes a substantial background to the A° ^ p 7U decays. 
However, the K*^^ will be used to cross-check the reconstruction algorithms and to clean 
the (and A^) samples. This is needed in order to avoid ambiguities when searching for 
the Aj, J/\|/ A° decay (and charge conjugate). 

6.2.1 The K*^s^ ^ signal in the region near the target/SSDs 

For decays in this region, the two pion SSD-PWC linked tracks are required to 
form a vertex (distance of closest approach < 50 |J,m). In order to resolve the vertex 
from the primary vertex the K^^ decay vertex should be at least 0.2 cm downstream of the 
J/\|/ decay vertex (the assumed production point of the K^^). Again, according to Monte- 
Carlo, an impact parameter of a track to the primary vertex of less than 200 |i.m means the 
track is associated with the primary vertex. Thus, the impact parameter of each pion 5^ 
relative to the K^^ production point must be greater than 200 |i.m. Also according to 
Monte-Carlo the impact parameter of the reconstructed K^^ relative to its production point 
must be less than 120 |i.m. A clear signal for the K*^^ is observed at the nominal K^^^ mass 
of 0.497 GeV/c, as seen in Figure 6.11. 
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Figure 6.11. Invariant-mass distribution using only SSD-PWC tracks for % % track 
combinations intercepting in the target/SSDs region. 
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6.2.2 The K^^— > n signal in the regions upstream and inside 
the dipole magnet 

In both of these regions the requirements are that the pion PWC tracks are not 
linked to any SSD track, and that the minimum relative pp is less than 0.4, as in the A*^ 
case. However, unlike in the — > p 7U decay where the proton tends to follow the 
direction of its parent A^, the pions from the K^^ decay have no prefered direction. Thus, 
no impact parameter is imposed on the pions. To enhance the k'^^ signal its transverse 
momentum is required to be greater than 0.4 GeV/c. These cuts resulted in clear signals in 
both regions for the K^^ at the nominal mass of 0.497 GeV/C. See Figure 6.12 and Figure 
6.13. 
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Figure 6.12. Invariant-mass distribution using only PWC tracks for 7t % combinations 
for tracks intercepting in the region upstream of the dipole magnet. 
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Figure 6.13. Invariant-mass distribution using only PWC tracks for 7t % combinations 
for tracks intercepting in the region inside the dipole magnet. 
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6.2.3 Distinguishing the A — > p 7i" (and charge conjugate) 
decays from the k"^ — > Ji* n decays 



The — > 7t^ 7t and A° — > p 7t (or — > p 7t^) decays are difficult to 
distinguished them from one another without particle identification. One can, however, try 
to distinguish the two decays by plotting the Podolanski-Armenteros plot [51], defined as 
the plot of the magnitude of the pj of either decay daughter particle relative to the direction 
of flight of the parent particle versus the asymmetry in of longitudinal momentum (in the 
laboratory frame) of the daughter particles, (pl - Pl)/(pl + Pl) • ^^^Y show that 
there is a kinematically allowed region bounded by: 



Pl-Pl E -E 
pt+PL E++E 

2p 
E+ + E 



where E , E , and p are the energy and momenta of the daughter particles in the 

* 

center-of-mass frame, and Pxmax - P • Figure 6.14 shows the Podolanski-Armenteros 
boundaries for the K^^ — > 7t^ 7t , A^ — > p 7t , and A^ — > p 7t^ decays, and Figure 6.15 
shows the Podolanski-Armenteros distribution of the A*^ (and A^) candidates decaying 

2 

upstream the dipole magnet. These include all the combinations between 1.105 GeV/c and 

2 

1.125 GeV/c in Figure 6.7. A cut on the pj of the decaying particles relative to the 
direction of flight of the parent particle to be less then 0.18 GeV/c is imposed to reject 
possible -^Ti^ % events. This value was suggested from Monte-Carlo since 100 % of 
the A^ ^ p 7U , and A^ — > p decays are retained. The A^ (and A^) candidates were 
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reduced to 10, 471, and 578 events, in the target/SSDs, upstream of the dipole magnet, and 
inside the dipole magnet regions, respectively. 
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Figure 6.14. Podolanski-Armenteros plot for the k!^^ — > tt^ tt , A*^ ^ p ti , and A^ ^ p 
decays. 
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Figure 6.15. Podolanski-Armenteros plot for the A*^ (and A^) candidates decaying 

2 

upstream the dipole magnet. These included all the combinations between 1.105 GeV/c 

2 

and 1.125 GeV/c in Figure 6.7. The lines are the kinematical boundaries for the K ^ — > 
tC Ti , ^ p 7U , and A^ — > p decays (see previous Figure). 
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However, even after applying the cut described in the previous paragraph there is 
still — > 7t^ 7U contamination from the outer sides of the and boundary regions. 
Figure 6.16 shows the mass spectrum of p 7t and p7t^ combinations that passed the 
Podolanski-Armenteros cut from the region upstream of the magnet, when both tracks are 
assigned pion masses. The solid Une is the expected distribution obtained from Monte- 
Carlo, normalized to have the same number of events as the data in the in the mass interval 

2 2 

between 0.3 GeV/c and 0.45 GeV/c . A similar plot is shown in Figure 6.17 for those 
events in the region inside the magnet. 

There is a clear K^^^ ii' n contamination. Because of this K^^, n 
background, any A° (or A^ ) candidate that has an invariant-mass consistent with that of the 
under the assumption of both tracks being pions, is excluded from the A° (and A^) 

2 

sample. The K j, mass range in each of the regions is: (0.48 - 0.51) GeV/c , (0.475 - 
0.525) GeV/c^, and (0.45 - 0.55 GeV/c^), for the target/SSD, upstream of the magnet, and 
inside the magnet, respectively. This cut retained 90% (9 combinations), 67% (318 
combinations), and 49% (285 combinations), in each region respectively. The efficiency 
of this cut according to Monte-Carlo is 90%, 67%, and 49%, for each region, respectively. 
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Figure 6.16. Invariant-mass spectrum of p 7U and pn combinations that passed the 
criteria of being a decaying in the region upstream of dipole magnet when both tracks 
are assigned pion masses. The solid line is the expected distribution obtained from Monte- 
Carlo and the dots are from data. 
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Figure 6.17. Invariant-mass spectrum of p tt and pn combinations that passed the 
criteria of being a A*^ decaying in the region inside of dipole magnet when both tracks are 
assigned pion masses. The solid line is the expected distribution obtained from Monte- 
Carlo and the dots are from data. 



CHAPTER 7 



THE J/\|/ A" DECAY CHANNEL 

This chapter includes a description of the J/\|/ A*^ Monte-Carlo simulation, as 
well as the search method and the results from the search for A^— > J/\|/ A^ events in our 
experiment. 

7.1 The Ab^ J/\|/ a" Monte-Carlo simulation 

A Monte-Carlo simulation of A,,— > J/\|/ A*^ events was used to tune the 
reconstruction program to search for such events in our data. It was also used to measure 
the acceptance of our detector and the efficiency of the reconstruction program for this t5^e 
of events. The Monte-Carlo was a GEANT3-based full detector simulation [52]. 

The Monte-Carlo simulation generates a bb pair according to the next-to-leading 
order (NLO) calculations of Mango, Nason, and Rodolfi (MNR), and it includes the 
following parameters: the mass factorization scale Q = m^J2, the Aq^q = A5 = 204 MeV, 
and the MRS235 and SMRS parton distribution fuctions for the nucleon and pion, 
respectively [21]. The two-dimensional distribution of Xp vs. p^ for the generated b 
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quarks is shown in Figure 7.1, and the projections on Xp and are shown in Figure 7.2. 
The s are assumed to have the same momentum as their parent b quarks. They are 

-12 2 

assigned a mean lifetime of 1.07 10 s, and a mass of 5.641 GeV/c . Both of these 
values are taken from the PDG book [12]. In each event, one of the A^, s in chosen at 
random and forced to decay into J/\|/ A^, with J/\|/ — > and A^ — > p 7U , or the charge 
conjugate reactions in the case of a A^. It is assumed that the A^, s are unpolaraized and 
that their decay is isotropic in the center-of-mass frame. The other b hadron was forced to 
decay into a kaon plus a randomly chosen number of pions ( up to a maximum of 5). With 
the remainder of the energy, hadrons in the underlying event are generated subject to 
energy, momentum and charge conservation, according to a longitudinal phase space. See 
reference [53] for a detailed description of the model. 
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Figure 7.2. (a) Generated Xp distribution of the b quark from Monte-Carlo events 
(projection from the two-dimesional distribution); and (b) generated distribution of the 
b quark from Monte-Carlo events (projection from the two-dimesional distribution). 
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All the events in the generated file were run through the detector geometry by 
GEANT. The resultant two-dimensional dimuon acceptance as a function of Xp and 
from A|,— > J/\|/ A^ events is shown in Figure 7.3. Our detector only accepted dimuons in 
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the Xp> region. The integrated acceptance for these J/\|/ s is 0.18 in the region Xp > 0. 

The events that had a dimuon accepted by the apparatus were digitized as 
spectrometer hits with appropriate detector noise and efficiency to produce a fake raw-data 
file. This file was then run through the same reconstruction and analysis program as that 
used for the data. Figure 7.4 shows a comparison of the track multiplicity in the Monte- 
Carlo with that in the data. Figure 7.5 and Figure 7.6 compare the momenta and p^ of 
hadron tracks from Monte-Carlo and data. 




Figure 7.4. Monte-Carlo and data distributions of the reconstructed track multiplicity. 
The solid line represents the Monte-Carlo and the dots the data. 
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Figure 7.5. Monte-Carlo and data distributions of the momenta of charged hadron tracks. 
The soUd line represents the Monte-Carlo and the dots the data. 




0.5 1 1.5 2 2.5 3 3.5 4 

(GeV/c) 

Figure 7.6. Monte-Carlo and data distributions of the of charged hadron tracks. The 
soUd line represents the Monte-Carlo and the dots the data. 
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Figure 7.7 shows the reconstructed background subtracted signals for the J/V|/ from 
Monte-Carlo and data. Figure 7.8, Figure 7.9, and Figure 7.10 show the reconstructed 
background subtracted signals of the (and A^) in the regions near the target/SSDs, 
upstream of the dipole magnet, and inside the dipole magnet, respectively. As seen, the 
resolution of the reconstructed J/\|/ mass and A° (and A^) from Monte-Carlo are in good 
agreement with those from the data. 
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Figure 7.7. The reconstructed J/\|/ signal after background subtraction, the solid line is 
from Monte-Carlo and the dots are from data. The Monte-Carlo is normalized to have the 
same number of events in the signal region as the data. 
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Figure 7.8. The reconstructed A° signal after background subtraction for those A*^ s that 
decayed in the target/SSD region. The solid line is from Monte-Carlo and the dots are from 
data. The Monte-Carlo is normalized to have the same number of events in the signal 
region as the data. 
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Figure 7.9. The reconstructed signal after background subtraction for those A^s that 
decayed upstream of the dipole magnet. The solid line is from Monte-Carlo and the dots 
are from data. The Monte-Carlo is normalized to have the same number of events in the 
signal region as the data. 
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Figure 7.10. The reconstructed signal after background subtraction for those A^s that 
decayed inside the dipole magnet. The solid line is from Monte-Carlo and the dots are from 
data. The Monte-Carlo is normalized to have the same number of events in the signal 
region as the data. 



106 

7.2 The J/\|/ a" search method and results 

The topology of the J/\|/ A° decay is illustrated in Figure 7.11. 

Combinatorial background to the A,, signal comes from three sources. The first one is 
from p 7C~ and \x^\x~ backgrounds underneath the A^ and J/\|/ signal regions, respectively, 
(see Figure 5.1, Figure 6.3, Figure 6.7, and Figure 6.10). The second is from real J/\|/ s 
combined with real A*^ s in the same event to give a fake A^^. The third is by wrongly 
assigning a proton mass to one of the pion tracks so that a K^^ ii' n could look like a 
A*^ ^ p 71 (or A^ — > p 7C^), since the experiment did not have particle identification 
capabilities. 



To search for the A^ J/v|/ A decay, the (9 + 318 + 285 ) A (and A ) candidates 
that survived the K^^ mass cut (see section 6.2.3) were required to have their transverse 

Dipole magnet PWCs 




Figure 7.11. The topology of the J/\|/ A decay. 
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momenta greater than 0.8 GeV/c. This cut was imposed to filter out A s produced in the 
primary interactions. Since A^s originated from A|,s have an average transverse momenta 
of 1.3 GeV/c whereas the directly produced A^s had on average transverse momenta of 0.6 
GeV/c (see Figure 7.12). This cut reduced the data to 2, 86, and 61 a'^ (and A^) 
candidates, in the target/SSDs, upstream of the dipole magnet, and inside the dipole magnet 
regions, respectively. According to Monte-Carlo this cut was 67 %, 75 % and 85 % 
efficient in the three regions, respectively. 

The A^ (and A^) candidates that survived the transverse monenta cut were then 
checked to see if a J/v|/ candidate that passed the dimuon refit (see chapter 5) and did not 
have its vertex upstream of the primary vertex, was present in the event. This criteria 
reduced the A^ (and A^) data sample to (0 + 7 + 8). Table IX shows a summary of the 
number of A*^ (and A^) candidates that survived each cut employed to reduce to this 
number of candidates. 



Table IX Summary of the number of A (and A ) candidates that survived each cut 
employed to select the sample used to search for the Aj, — > J/\|/ A° decay. 



Region 


Initial 
sample 


Podolanski-Armenteros 
Cut 


K 5 mass 
Cut 


A -Pt 
Cut 


J/\|/ 
requirement 


I 


10 


10 


9 


2 





II 


479 


471 


318 


86 


7 


III 


582 


578 


285 


61 


8 
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Figure 7.12. (a) A -p^ distribution from data, and (b) from the Aj,— > J/\|/ A Monte-Carlo. 
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The invariant-mass for the (0 + 7 + 8) J/V|/ A combinations summed over the three 
reconstruction regions is shown in Figure 7.13 (a), assuming the cuts described above. 
The distribution of invariant-mass of A^, s reconstructed from Monte-Carlo simulation is 

2 

shown in Figure 7.13 (b). This shows the A,, signal region to be in the range 5.4 GeV/c to 

2 

5.9 GeV/c . It is easy to see from Figure 7.13 (a) that there are two events in this mass 
range. Thus, there are two A^, — > J/\|/ A^ candidates. To estimate the shape of the 
combinatorial background to this distribution, a "wrong-frame" J/V|/ A^ background was 
produced by combining real J/\|/ s with real A° s (and A^ s) from different events and 
applying the same set of cuts as for the signal. The distribution of invariant mass from 
"wrong-frame" events is shown in Figure 7.13 (c). It is clear that the signal region is not 
in the same region where the wrong-frame background peaks. 

In principle with the data shown so far one can also ask about the b'^ ^ J/\|/ K° 
decay. However, due to the smaller Br(B*^ J/\|/ K°) (= (7.5 ± 2.1) x 10""^ [12]), 
compared to Br(Ay, J/\|/ A°) (= 1.8 x 10"^ [7]), and F(b ^ K° ) (» 0.2) compared to F(b 
— > Aj, ) (~ 0.1), and considering we can only observe K^^ and not k\, the expected 
number of B° — > J/\|/ events is about 1/24 of those expected from Aj, — > J/\|/ A° 
decays^ Thus, we do not expect to see this decay in our experiment. 



' See chapter 8 for the expected number of A}, — > J/\|/ A decays. 
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Figure 7.13. (a) The invariant mass for the J/\|/ A combinations; (b) the invariant mass 
for the J/\|/ combinations from Monte-Carlo; (c) J/\|/ A^ wrong frame background. 



CHAPTER 8 



EVALUATION OF THE A^^ J/\|/ A" PRODUCTION RATE 

This chapter gives a detailed explanation of the detection efficiencies for the J/\|/s 
and A*^ s, and describes the procedure used to determine an upper limit to the Aj,— > J/\|/ A*^ 
production rate, F(A|,) * Br(A(, J/\|/ A°). The upper limit is calculated based on two 
events found in the A^, region (see section 7.2). It is then recalculated twice using a 
modified set of cuts. 

8.1 Acceptances and efficiencies 

The acceptances and efficiencies to reconstruct each signal were measured using the 
Monte-Carlo summation described in section 7.1. The acceptance for J/V|/ — > |i with Xp 
> from Aj, events is 18 % (see section 7.1). The reconstruction efficiency is 59 % and it 
breaks down in the following way: 

i. The efficiency of linking each muon throughout the entire detector is 86 %. 



Ill 
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ii. The invariant mass cut of 2.85 GeV/c^ < M^^ < 3.35 GeV/c^ retains 93% of the 
J/\|/s. 

iii. The dimuon refit to the J/\|/ explained in chapter 5 and appendix A is 92 % efficient. 

iv. The requirement that the be vertex is not upstream of the primary vertex is 93 % 
efficient. 

2 

Thus, the reconstruction efficiency for J/\|/ — > it |X" from events is (0.86) 
(0.93) (0.92) (0.93) = 0.59. The "PSl" dimuon preselection described in section 4.1.1 is 
71 % efficient for J/\|/s from A|, s. Therefore, the total detection efficiency for J/\|/ s from 
AfeSis e^^VfromA^, ^ (0.18)(0.59)(0.71) = 0.075. 

The kinematical acceptance, that is the total fraction of the A^s that decay in the 
three regions once the J/\|/ has been accepted is 63%. The individual kinematical acceptance 
for each region is 3 %, 25 %, and 35 %, for the target/SSD, upstream of the dipole 
magnet, and inside the dipole magnet regions, respectively. The geometrical acceptance 
times reconstruction efficiency of A^ s in each region is 8 %, 16%, and 29 %, respectively. 
After including the K^^^ mass cut which is 90 %, 67% and 49% efficient ( see section 
6.2.3), and the A^-p^ cut on A^ (or A^) which is 67 %, 75 % and 85 % efficient (see 
section 7.2), the geometrical acceptance times reconstruction efficiency of A^ s in each 
region reduces to 5 %, 8 %, and 12 %, respectively. Therefore, the weighted average total 
acceptance times the reconstruction efficiency for the A° once the J/\|/ has been accepted is 
8^° = (0.03) (0.05) + (0.25) (0.08) + (0.35) (0.12) = 0.064. 
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8.2 Upper limit on the production rate 



The observed number of events in the Aj, signal, N^jg, is related to the exclusive 



branching ratio Br( Aj, — > J/V|/ A*^ ) through 



where vfromAb ^.j^^ ^^^^j detection efficiency for J/\|/ s originating from A^, s; 8^ 
is the total A^ detection efficiency; F(A^,) is the production fraction (i.e., the probability of a 



b quark to hadronize into a or into another beauty baryon decaying to Aj,; Br(A(, — > J/\|/ 

bb 



A° ) is the branching ratio for — > J/\|/ A^; o, r- is the bb production cross-section; 



Br(J/\|/ |i+ |i ) is the branching ratio for J/\|/ ^ |i+ |i ; Br(A*^ ^ p;r") is the branching ratio 
for A*^ — > p7t and L is the integrated luminosity. The prompt J/\|/ signal from our beryllium 
target is used to compute L, this way there is a reduction of systematic uncertainties arising 
from efficiency corrections due to dead-time, muon halo, pretrigger and the dimuon trigger 
processor. L is written as: 



L = ^. ^ — Y, (8.2) 



8J^^aj/yBr(j/i|/^rt")^ 



where 8^^^ is the total detection efficiency for prompt J/\|/ s. In reference [11] we report, 
Nj/v,, = 9,800 ± 130 on Be with Xp> 0.1. This corresponds to Oj/^ Br(J/\|/ ^ |i )/A = 
9.2 ±1.2 nb/nucleon for Xp> 0.1 on Be [11]. The acceptance for these J/V|/ s with Xp > 0. 1 
is 43%, the reconstruction efficiency is 64 % [11], and the "PSI" dimuon preselection 
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described in section 4.1.1 is 79 %. Thus, the total detection efficiency for prompt J/\|/s is 
e^^^= (0.43)(0.64)(0.79) = 0.217. 

The total detection efficiencies gJ^VfromAb^ ^ gJ/v listed in Table X. 
The Br(J/\|/ ji") = 5.97 %, and Br (A^ ^ p 7t") = 63.9 % [12]. Also E672/E706 
measured a^:^ to be 47 ± 19 (stat) ± 14 (sys) nb/nucleon for Xp > 0.0 [11]. From 
equations 8.1 and 8.2: 

Nsig e^^^cJj/v Br(J/\|/ ^ ^i^ ^i") 

F(Ab)*Br(Ab^J/\|/A°)= . (8.3) 

Nj/^2 e^/'VfmmA^ gA" ^ gj. ^^0 ^ p ^-^ 

To compute an upper limit to F(A^) * Br(Aj, — > J/\|/ A^), the 2 events in the A^, 
signal region (see section 7.2) are treated as signal with zero background. According to 
Poisson statistics the maximum number of signal events that the 2 events can statistically 
fluctuate up to is 5.3 events at the 90 % C.L'. Therefore 

F(Ab) * Br(Ai, ^ J/\|/ A*^) < 6.2 x 10"^ at 90 % C.L. 

This does not contradict the value measured by UAl of F(Aj,) * Br(A|, — > J/\|/ A^) = 1.8 ± 
1.1 X 10 [7]. Using the UAl measurement and the set of cuts described in the previous 
chapters we expect to find 0.2 A^, — > J/\|/ A*^ events. To better understand and measure the 
F(A|,) * Br(Aj, — > J/\|/ A°) upper limit, this number was recalculated twice under a modified 
set of cuts than those used above. 

' See ref. [12] page 1279 



115 



Table X. The total detection efficiencies for the different particles, which are the products 
the kinematical and geometrical acceptance, reconstruction efficiencie and preselection 
efficiencie. 





Acceptance 
(%) 


Reconstruction 
efficiency (%) 


Preselection 
efficiency (%) 


Total detection 
efficiency (%) 


J/\|/ s from A(,s 


18 


59 


71 


gJ/\|/fromAb _ 7 5 


prompt J/\|/s 


43 


64 


79 


8^'^ =21.7 


A°(and) 




6.41 




= 6.4 



The first set was made without including the K ^ mass cut mentioned in section 
6.2.3. Without this cut the A° total detection efficiency 8^^" = (0.03) (0.05) + (0.25) 
(0.12) + (0.35) (0.25) = 0.119 (see section 8.1 for the definition of each factor). 
Consequently removing this cut, the A^ total detection efficiency increases approximately 
by a factor of 2 and therefore increases the expected number of A;, signal events to about 
0.4 events, maximizing the signal. Without the k'^^ mass cut, the A^ (or A^) data sample 
contains 10, 471, and 578 candidates in the target/SSD, upstream of the dipole magnet, and 
inside the dipole magnet regions, respectively (see section 6.2.3 for all other cuts applied). 
Applying the cut A*^-pj to be greater than 0.8 GeV/c reduces the data sample to 2, 122, and 
162 candidates in the three regions, respectively (see section 7.2 for the motivation of this 
cut). After searching each event in this data sample for a J/\|/ candidate that passed the 
dimuon refit (see chapter 5) and did not have its vertex upstream of the primary vertex, the 
A° (and A^) data sample reduced to (1 + 22 + 25) events. The invariant-mass for these 48 
J/\|/ A^ combinations is shown in Figure 8.1 with the solid line. The dashed line represents 

' This value is the weighted average acceptance times reconstruction efficiency for the 
A (and A^) in the three reconstruction regions combined, once the J/\|/ has been accepted. 
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wrong-frame background normalized to have the same number of events in the region 

2 2 

between 6.5 GeV/c and 9.0 GeV/c as the data. Recall that the A}, signal region is the 

2 2 

mass interval between 5.4 GeV/c and 5.9 GeV/c . There are 7 A^, candidates in this signal 
region (number of events observed), and 12 background events in the wrong-frame 
background (number of expected background events). According to Poisson statistics, the 
maximum number of signal events that 7 observed events with an expected background of 
12 events, can fluctuate up to is 4.2 events at the 90 % C.L^ Using equation 8.3 were 
every value is the same expect for 8 which is now 0. 1 19, one finds that 

F(Ab) Br(Ab J/V|/ a") < 3.1 x 10"^ at 90 % C.L. 

The second set of cuts also does not included the K^^ mass cut, but in addition the 
J/V|/s are required to originate from a vertex downstream of the primary vertex in the event. 
In this set of cuts the background is reduced to a minimum. The reconstruction of J/\|/ s 
with xp > from A^ events is now (0.86)^ (0.93) (0.92) (0.60) = 0.38 (see section 8.1 for 
the definition of each factor). Thus, e^^^^"""^" = (0.1 8)(0.38)(0.71) = 0.048. Giving an 
expected number of event of about 0.2. Requiring that the J/\|/ s emerge from secondary 
vertices, reduces J/\|/ A° data sample to (1 + 1 + 2) events. The requirements for a J/\|/ to 
be originated from a secondary vertex are the same as those used in the J/V|/ selection 
process for our bb cross-section calculation in [1 1]. (i) The primary vertex in the event is 
required to have at least 3 SSD-PWC linked tracks associated with the vertex, (ii) The J/\|/ 
vertex is required to be at least 2.5 mm downstream of the primary vertex, and (iii) a 



' See ref. [12] page 1279 
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significance greater than 3 is required for both longitudinal and transverse separations 
between the primary and secondary vertex, with the significance defined as the separation 
divided by the combined uncertainty. The invariant-mass for the (1 + 1 + 2) J/\|/ 
combinations is shown in Figure 8.2. There are no events in the Aj, signal region. For the 
purpose of computing an upper limit to F(A|,) * Br(Ai, — > J/\|/ A°), the signal is zero and 
the is zero background. According to Poisson statistics, the maximum number of signal 
events that zero signal and zero background can statistically fluctuate up to is 2.3 events at 
the 90 % C.L\ Using equation 8.3 as before but with e-'^^^^^^t = 0.038, one gets 

F(Ab) * Br(Ab ^ J/V|/ A°) < 3.2 x 10"^ at 90 % C.L. 



' See ref. [12] page 1279 
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Figure 8.1. The invariant mass distribution (soUd) for the J/\\f A combinations, without 
using the mass cut, and the wrong-frame background (dashed). 
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Figure 8.2. The invariant-mass distribution for the J/\|/ A° combinations, without using 
the mass cut, and requiring that the J/\|/ s originate from secondary vertices. 
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CHAPTER 9 



CONCLUSIONS 

An extensive analysis was performed on 8.0 pb %f dimuon data produced in 7t A 
collisions at 5 15 GeV/c to search for Aj, events in the decay channel — > J/\|/ A°, with J/\|/ 
— > |X and A° — > p 7U (and for the conjugate reactions). A refit to the muon tracks from 
J/\|/ decays was performed with the mass constraint of two-body decay and with the 
constraint that both muon tracks intersect at a common point. The A*^ s were identified by 
their characteristic decay, giving the larger fraction of their momenta to the protons. Using 
an iterating algorithm, A° s were reconstructed in three regions of the E672/E706 
spectrometer, near the target/SSDs, between the SSDs and the dipole magnet, and inside 
the dipole magnet. A total of 575 ± 35 A^ (and A^) background subtracted events were 
reconstructed. A K^^ — > 7t^ 7t signal was also reconstructed in all three regions using the 
same technique. The K^^ mass signal was used to cross-check the reconstruction algorithm 
and to clean the A^ (and A^) data sample, giving (9 + 318 + 285) A° (and A^) candidates. 

To search for the Aj, — > J/\|/ A*^ (and charge conjugate reaction), J/\|/ s that passed 
the muon refit were combined with clean A^ s (or A^ s) when they existed in the same 
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event, giving a total of (0 + 7 + 8) J/V|/ A (and charge conjugate) combinations. The 
results show 2 events in the Aj, mass region. Considering the two A}, events as signal with 
zero background, an upper limit to F(A^) * Br(A^ — > J/\|/ A^ ) was found to be less than 

-2 

6.2 xlO at 90 % C.L. An upper limit was also calculated, without using the K ^ mass cut 
for the A° s (and A^s), then, also requiring that the J/\|/ s originate from secondary 
vertices, giving that F(A^) * Br(Ab J/\|/ A° ) < 3.1 xlO"^ at 90 % C.L and F(Ab) * 

-2 

Br(Ai, ^ J/\|/ A ) < 3.2 xlO at 90 % C.L, respectively. The upper limits should be 

-3 

compared to (1.8 ± 1.1) x 10 measured by UAl [7]. 

This study was performed using a Monte-Carlo simulation for bb production from 
next-to-leading order (NLO) calculations of Mango, Nason, and Rodolfi (MNR), which 
included the mass factorization scale Q = the Aqqj = A5 = 204 MeV, and the 

MRS235 and SMRS parton distribution functions for the nucleon and pion, respectively 
[21]. In addition It was assumed that the A|, s are unpolaraized, and no fragmentation was 
included. 



In conclusion E672 did not have the sufficient sensitivity to contradict the UAl 

.3 

measurement, nor give a lower limit than the ones established by ALEPH (< 0.4 x 10 at 
90 % C.L.), OPAL (< 1.5 xl0"^90 % C.L.), and CDF (< 0.5 x 10"^ at 90 % C.L.), (see 
Refs. [8] and [9]). However, ALEPH and OPAL both used e^e interactions, UAl and 
CDF usedpp interactions , and CDF had a limited sensitivity to reconstruct A^ s (and A^s) 
that have a n (and n') with pj < 0.4 GeV/c [9]. Here we used n A collisions and 
reconstructed A° s (and A^s) with the p^ of the K (and n^) between GeV/c and 0.4 
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GeV/c. Thus, our search for the Aj, J/v|/ A (and charge conjugate reaction) in 
complementary to that of the others mentioned above. 



APPENDIX A 



FIT TO THE J/\|/ 



The muons that forai the decaying J/\|/ are fitted using a Least-Squares method 
(explained in detail below) with the constraint that the two muons intersect at a common 
point in space and that the invariant mass of the combined pair is equal to 3.097 GeV/c ^ . 

2 

The X of the fit is used to reduce the combinatorial background when there are more than 
two muon candidates in a given event. The fit improves the resolution of the measured 
momenta of the muons by a factor of 2 and vertex position by 15%. 



General least-squares estimation 

The mathematical formulation of the 
making any reference to any special physical 



with constraints 

iterative procedure will be derived without 
problem. 



' This is the value published by [12]. 
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The iteration procedure 

Let M = {mj,m2,...,mj^}be a vector of N measurable variables, which have initial 
measurements M° = |m|',m2,...,m° }, with errors contained in the co variance matrix 
E(M°). In addition, let U = {ui,U2,...,Uj} be a set of J unmeasurable variables. The 
M variables and the U variables are related and have to satisfy a set of K constraint 
equations 

f^(mi,m2,...,mj,,UpU2,.--jUj) = , k = l,2,...,K . 

According to the Least-Squares Principle, the best estimates of the M and U 
variables are those for which 

X^(m,M°) = (m''-m)^E"'(m°)(m°-m) = minimum, and (A.l) 

f(M,u) = 0, 

where the superscript T indicates the transpose of the matrix. 

The problem will now be solved using the Lagrange multiplier method. 
Introducing the K component vector A, = {A,j,A,2,...,A,k^} of Lagrangian multipliers, the 
problem can be rephrased by requiring 

X' (m, U, i ) = (m'' - m)^ E"' (m° ) (m° - m) + 2 f (m, u) = min imum . ( A.2) 
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There are now a total N+J+K variables. The values of M , U and X that make the 

2 

X minimum, must also satisfy the following set of equations 

Vm%'=-2E"'(m°)(m°-m) + 2F;^ i = 0, (N equations) 

Vu x' = 2 i = , (J equations) (A.3) 

V^X' = 2 f (m, u) = 6 , (K equations) 

where the matrices Fm,Fu of dimensions KxN and KxJ respectively are defined by 



(F ) =^ (F ) =^ 



Simplifying equations (A.3) gives 

-E"' (m° )(m° - m) + F^i; = 6 , ( A.4a) 



Fix = 6, (A.4b) 
f (m,u) = 6 . (A.4c) 



The solution to the set of equations (A.4) for the N+J+K variables must be, in 
general, found by iterations producing successively better approximations. To do this, 
suppose that the n-th iteration has been performed and it is still necessary to find a better 
solution. For the n-th iteration the approximate solution is given by the values of 
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M",U",A,", corresponding to the function value (X^)". By performing a Taylor expansion 
of the constraint equations (A.4c) about the point (M°,U°) one obtains: 



af. 



df 



fk + L ^ (mr'-m:)-£ ^ (uf^-u;)+...= 0, k = l,2,...,K 



Neglecting all terms of second or higher order, this equation can be written as 



r + f;:j(m"^' -m")+f^(u°^' - u") = 6 , 



(A.5) 



where all superscripts n indicate that the f",F5^,Fu are to be evaluated at the point 
(M",U"). The other two equations in A.4, a and b, are then, 



E-'(M°+i-M0) + (F;^)"r"'=6, 



(A.6a) 
(A.6b) 



These two equations, together with the expanded constraint equations (A.5) make it 
possible to express all variables of the (n+l)-th iteration in terms of quantities of the 
preceding iteration. These solutions are 



U»+i = U" _ (pT 3-1 j-i pT 3-1 r ^ 

r^i = s-'[F+Fu(u"^^-u")], 



(A.7) 



T in+l 



M"^^ = M"-EFi^ I' 
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where 

r = r+F;;i(M°-M°), 
s-f;e(f^)\ 

In equations (A.7) the matrices F^,Fu,S, and the vector f are evaluated at the point 

J- J- Q 

(M",U"). The value of U can be found from the set of the constraint equation, setting 
f(M°,U°) = and solving for U°. With the new values for M"^',U"^'and T^' one can 
calculate the value of the function (x^)"^^ for the (n+l)-th iteration and compare it with the 
previous value of (x^)", by using 

(^'Y^' = (x-')' s r^'+2 (r^'f (a.8) 

where the matrix S is evaluated for the n-iteration. Once a satisfactory solution is found the 

2 

iterative procedure is stopped. This is when the value of the >C and the values of the 
M and U vectors are converging [54]. 

Kinematical analysis of the J/\|/ decay 

The following section applies the previous formulation to the J/\|/ — > [L^[L decays 
as seen in the MWEST spectrometer. Let us assume that the two muon tracks have been 
measured and that the track reconstruction program has provided for each muon track a 
first approximation for the kinematical variables: 1/p ( inverse momentum); m^ (track slope 
in the x-z plane ); my (track slope in the y-z plane ); b^ (intercept of the track with the x-axis 
in the x-z plane); and by (intercept of the track with the y-axis in the y-z plane); as well as a 
covariance matrix (E) for these values. Since the decay vertex of the J/\|/ (x, y, z) is 
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unspecified and the magnitude and direction of its momentum is unknown, the problem in 
this case involves six immeasurable variables 



U = j^,m^j/^,myj/^,x,y,z| 



and ten measurable variables for both muon tracks. 



r 

M = 



^— ,m ,,m ,,b ,,b ,m _,m _,b _,b _ 

„ yn+' XH+' yH „ xn ' yn ' x^ ' yH 



The algebraic constraint equations are: the three equations describing momentum 
conservation; the equation from energy conservation; and the geometrical line-equations of 
the tracks in the x-z and y-z planes. This gives a total of eight equations: 



m 



x J/\|f 



+ ■ 



m , 



m 



■ + ■ 



x|i 



Pj/\l/ 



^ m^ , + m^ , + 1 ^^Jm^ _ +m^ _ +1 ' 

xn yn „ Y xji yn 



f2 = 



m 



yJ/v 



■ + ■ 



m + 
yn 



m 



■ + ■ 



yti 



Pj/v 



mx,/v + myj/v+l 



m , + m , + 1 

xn yn 



m" +m" +1 ' 

xn yn 



f3 = 



■ + ■ 



■ + ■ 



m 



2 I 2 I 1 



P.^ 



m , + m , + 1 

x|X xn 



m + m +1 

X|i X\L 



f4 = 



1 2 

' 1 ^ 



1 2 

7 — V^'^l + 
1 



vPj/vy 

f , = -X + z m + + b + 



1 2 
P ^ 
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f, =-x + zm _+b 



XH- 



f7 = -y + zm ++b + 



f„=-y + zm _+b 



where nij/^ and are the masses of the J/\|/ and muon, respectively. Since the problem 
involves 8 constraint equations and 6 unmeasurable variables this is then a 2c-fit. 



Results from the fit 

In a given event, all the unlike-sign muon pairs with fully linked SSD-PWC- 

2 

MUON tracks, and with an invariant-mass in the mass interval between 2.85 GeV/c to 

2 

3.35 GeV/c , and having a distance of closest approach between the two muon tracks less 

2 

then 50 |J,m, were fitted with the above technique. For convergence, the values of the x , 
M, and U variables were required to be within 0.001 of the preceding values of each 

2 

variable, respectively. In each event only the dimuon with the lowest % was kept. Figure 

2 

A. 1 shows the distribution of the % per degree of freedom of the fit. A cut on the 

2 

X /d.o.f less than 5.0 was established. Before the fit, there were a total of 13,053 muon 
pairs. After the fit 12,340 dimuons survived the fit requirements. These were identified as 
J/V|/ s. 

From the — > J/\|/ A^ Monte-Carlo, the reconstruction efficiency of the initial track 
finding for J/\\f s from Aj, s was found to be 49.8 %, and by using the refit procedure this 
efficiency was improved to 59 %. In Figure A. 2 the dashed line shows the residuals of 
the J/\|/ momentum between the generated value from — > J/\|/ A^ Monte-Carlo and that 
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reconstructed by the initial track finding. The solid line shows the residuals of the 
momentum after the refit. Figure A. 3 shows the residuals of the z-coordinate of the J/\|/ 
decay vertex between the generated value from Monte-Carlo and the reconstructed, the 
dashed line for the initial track finding and the solid line for the refit. As one can see, 
refitting the J/\|/ resulted in an improvement of its momentum measurement resolution by a 
factor of 2, and in the vertex resolution by 15%. 



250 
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J/v|/ momentum residuals (GeV/c) 

Figure A. 2. Residuals of the J/\|/ momentum between the generated momentum value in 
the A),^ J/\|/ A^ Monte-Carlo and the reconstructed value. The dashed line is for the 
values reconstructed by the initial track finding. The solid hne shows the residuals of the 
momentum after the refit. 
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Figure A. 3. Residuals of the z-coordinate of the decay vertex between the value generated 
in the Monte-Carlo and the value reconstructed. The dashed line is for the values found by 
the initial track finding and the solid is for those after the refit. 



APPENDIX B 



MOMENTUM MEASUREMENT IN THE REGION UPSTREAM 

OF THE DIPOLE MAGNET 

Assuming the decay vertex of the A° (x, y, z) and the parameters of the PWC track 
^xD (x-slope), (y-slope), (x-intercept) by (y-intercept) are known, and that the 

J- I -Lp 

magnetic field is uniform throughout the dipole magnet, B = - B j . The momentum vector 
is estimated in the following way: First, the coordinates of the PWC track at the center of 
the magnet are computed, 

Wer = 197.73 cm, 

'^center ~ ^center 

Then, to estimate the bend of the trajectory of particle (proton or pion), the x-slope 
m^u of the imaginary line in the x-z plane from (x, z) to (Xj.gjj(gp z^g^gj) is computed 

Y — Y 

m — ^center 
^ ^center 
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see Figure B. 1. M^u is the slope of the proton or pion trajectory at the A decay vertex. 
Since the magnetic field is along the negative y-axis (see section 2.5.2) the electric charge 
associated with the track is then. 



m 



xD 



m 



xU 



The magnitude of the momentum of the track in the x-z plane is: 





V 

Pt 

^kick 













where |P^^ | = 0.3|B|L, and |B|is the magnetic field (in Tesla), and L (in meters) is the 
length of the field along the z-axis. From P^^^ one can calculate the z component of the 
momentum: 



P. = 



1 + m 



xU 



and from this the x and y components: 



Px - Pz ^xU' 



and 



Py - Pz ^yD- 
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► Dipole Magnet 

J z 

Figure B. 1. The thin lens approximation is used to compute the momenta of the proton 
and pion tracks of those A° s that decayed upstream of the dipole magnet. 



APPENDIX C 



MOMENTUM MEASUREMENT IN THE REGION INSIDE THE 

DIPOLE MAGNET 

Assuming that the decay vertex of the (x,y,z), and the parameters of the PWC 
track m^p (x-slope), (y-slope), (x-intercept), by (y-intercept) are known, and that 
the magnetic field is uniform throughout the dipole magnet, B = - B j. The magnitude of 
the momentum is determined by measuring the radius R of the circular trajectory described 
by the proton or pion particle inside the magnetic field. R is related to the momentum of the 
particle in the x-z plane by, P^^ =0.3 B R, where B is measured in Tesla, and R in 
meters. The momentum vector is estimated in the following way: Labeling the center of the 
circular trajectory as (x^, z^), see Figure C. 1, it is easy to see that, 

{x-xJ+(z-zJ=R^ 

(Xend - ^of + (Zend " Zo)' = J ^nd 

. Xend - Xq . 
Zend Zq 

where, Zg^^ = 318.53 cm, and Xg^^j = mxD+ b^. 
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(x ,z ) 
^ o' 




PWC proton or pion track 



Tan" (nixD) 



Dipole Magnet 



► 

y z 



Figure C. 1. Circular trajectory described by of a proton or pion of a A that decays inside 
the dipole magnet. 

These three equations have three unknowns, x^, z^, and R. It is not so difficult to 
solve for these unknown variables and obtain 



~ ^end 



1 (x-x,„j'+(z-z^„j' 



2 (x-x,„d) 



m 



-(Z-Ze„d) 



xD 
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R 



(Zend-Zo)- 



Having this, m^u can now be computed. 



m 



xU 



z„ - z 



X - x„ 



and the electric charge of the particle is found. 



Q = 



FxD-mxu 



As mentioned previously the momentum of the particle in the x-z plane is P^^ = 0.3 B R, 
From one can calculate the z component of the momentum: 



1 + m 



xU 



and from this the x and y components: 



and 



Py - Pz l^yD- 
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